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ABSTRACT 

This thesis presents the results of an investigation of the lateral 
performance functions of an airplane by the analysis of the transient 
response of the airplane to finite pulse of control surface deflection. 

The method is an extension of pulse technique used by Clementson, 

Ref .(2), to measure the longitudinal performance functions of an airplane 
and further described by Seamane, Blasingmme, and Clementson, Ref.(3), for 
more general cases. In this thesis the successful extension to highly 
oscillatory lightly damped responses is shown. Also a possible technique 
of analyzing divergent responses ia indicated. 

The instrumentation and calibration neceseary to measure inputs of 
aileron and rudder deflection and outputs of time-rate of geometric yaw and 
roll are described in detail. The separate deflections of each aileron 
were added electrically to give the total differential aileron deflection 
thus obviating the necessity for adding these deflections during analysis. 
Microsyn signal generatore were used to measure aileron and rudder deflections. 
Rate gyros were used to measure the angular velocities in roll and yaw. 
These were the only outputs measured due to the limited time available. 

Determination of the performance functions from the transient response 
was accomplished by approximating the Fourier Integral of the time responses 
of the inputs and outputs useing the method of summing approximating triangles. 

The technique of measuring the transient response to finite pulse 
inputs is shown to be accurate, reliable and reproducible with a low level 
of uncertainties. | 

Special methods of analyzing lightly damped oscillatory modes are 
developed and compared. The recommended method will vary somewhat with the 


form of the response curve but generally speaking the most satisfactory 


“(i= 








method consists of approximating the Fourier Transform of the output over the 
early portion of the time response to a point beyond which the response is, 
within the accuracy of measurement, only a damped sinusoid. At this value 

of time the Fourier Tranaform of the oscilletory mode obtained analytically 
is added vectorially to the Fourier Transform of the initial portion of the 
response. 

Comprrison with computed performance functions based on calculated 
derivatives showe qualitative agreement but clearly shows that these 
computed performance functions may differ considerably from the actual 
performance functions of the full scale airplare in flight. 

The saving in instrumentation and flight time effected by the pulse 
technique over that required by the sinusoidal technique is a major 
factor in recommending the former. Changes in the mass and moment of 
inertia due to fuel consumption are negligible for the short duration of 


teet flights made possible by using the pulse technique. 
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OBJECT 


fo inetrument a B-chJ airplane and measure the 
lateral transient responses in rate of roll and rate of 
goometrio yaw to pulses of aileron and rudder deflection. 

To determine the lateral performance functions of 
the airplane from theee reeponses by approximating the 
Fourier Integrale of the input and output. 

fo evaluate this. method of obtaining these 


performance functions. 


-j)|- 





INTRODUCTION 

Recent demands for high-performance automatic control of aircraft 
for both military and commercial uses have focused attention on the 
airplane as the dynamic system being controlled. It is fundamental that 
the synthesis of a controlling system starts with the best possible 
evaluation of the dynamical system to be controlled. Additional interest 
in the dynamics of the airplane arises from the problems of stabilizing 
some of the auxiliary equipment utilized in modern airplanes, such as 
navigating and fire control equipment. 

The frequency response of an airplane can be calculated from computed 
derivatives based on wind tunnel tests. However, the comparison of these 
calculated responses with those measured on the actual airplane has shown 
appreciable discrepancies, due in part to the difference in stability 
derivatives obtained in the wind tunnel from those of theairplane, and in 
part to ignorance of the moments of inertia about the axes of motion. 

The most generally used methods of synthesis and design of automtic 
control systems requires the dynamic characteristics of the airplane to be 
presented as a frequency response. A measured frequency response is most 
often obtained vy the direct method of applying sinusoidal inputs and 
measuring the amplitude and phase angle of the outputs. This ia done for 
enough frequencies to define the complete frequency response. 

Pioneer work using the sinusoidal oecillation technique was carried 
out by the Cornell Aeronautical Laboratories with outstanding results. Pef:!-) 
These tests required an auto~pilot, a sinusoidal input generator and measur 
ing equipment. Extensive flight time, with its consequent expense, is 
associated with this technique. 

An alternative technique, that of measuring the transient response to 


a step function input, has the drawback that such an input to the elevators 


=2— 





or ailerons will in a very short time cause the airplane to rotate throgh 
a large angle and hence will have displaced it so far from the original 
position of equilibrium that the motion is nomlinear. Such motion at the 
present time can be analyzed only by extremely tedious and lengthy methods 
of numerical —— 

A technique which obviates a large part of the instrumentation and saves 
a great deal of flight time is that of measuring the transient response to 
a finite pulse of control surface deflection. This method was successfully 


(Ref .2.) in determining frequency response in the case of 


used by Clementson 
longitudinal motion. The mathematical background and the technique of 
application are presented in a paper by Seamans, Blasingame, and 
Clenenteon, 2°! *3+) 
In extending the pulse technique to the measurement of the lateral 
modes of motion of the airplane, the characteristic of spiral divergence, 
which exiatse in many airplanes, must be taken into account. The Fourier 


(Ref .3.) are valid only for 


Transformation method and its approximation 
stable systems. Therefore, it becomes necessary to identify and remove 
the unstable mode prior to transformining the pulse response. Furthermore, 
the oscillation that occurs in lateral motion is lightly damped. This 
oscillatory mode must be given special treatment, or the analysis will 
ol exceedingly tedious, in the regions of resonance. 

Methode of analyzing the oscillatory mode were successfully used in 
thia investigation to find the lateral performance functions of the B—25a 
airplane. Howevor, this airplane did not have any measurable characteristic 


of spiral divergence, therefore, to illustrate a method of analysis when 


euch a characteriatic is involved, a simple system was analyzed. 
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of spiral divergence, therefore, to illustrate a method of analysis when 


such a characteristic is involved, a simple system was analyzed. 





CHAPTER 1 


DETERMINATION OF AIRPLANE LATERAL MOTION PERFORMANCE FUNCTIONS BY PULSE 
TECHNIQUES 


The airplane in flight may be considered as an operating component 
which operatee on epecified inputs to produce certain outputs. More than 
oné output may result from a given input. The operation on an input to 
produce a specified output is expressed mathematically by a concept of 
wide generality called a Performance Operator. 

By restricting the airplane and its controls to small excursions 
from an equilibrium condition, its motion may be defined by a system of 
linear differential equations. Then the reeponse can be related to the 
input by means of a Performance Function. The concept of Performance 


Operator and Performance Function were both developed by Dr. C. &. Draper 


and his associates of the Instrumentation Laboratory of the Maseachusetts 
Inetitute of Technology. A complete discussion of their application to 
airplane dynamics is contained in Refs. (1) and (10). The relation between 
the two may be ‘atkarion by the following equation: 


+jJwt _  +jwt 


When the output of an operating component can be related to the 
input by a set of linear differential equatione, and the input ie 
described by a complex exponential, the output is described by 

jt, t + PA) 
Tout) * Wouta © 
The ratio of the output to the input then defines the performande 


function of the operating component relating the two. In equation form 
jLwt PA) oat) 


q * 
(PF) ¢ Man Vout? . itv, t + (Pa), ] 


Tina 5 
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jw it 
The term e f factors out, and it is seen that the performance 


function is independent of time. For use in the synthesis of control 
or stabilising systems, the performance function is usually presented 
graphically, either as a locus on the complex plane or as an amplitude 
ratio and phase angle plotted versus forcing angular frequency on the 
real plane. (1)(2)(3)(10) 

It is the need for the frequency spectrum of the airplane’s dynamics 
that has extended the study of theese dynamics beyond the classical approach 
which, until the last few years, has been concerned with only the 
characteristics equation of the dynamic response and sensitivities. This 
approach does not give sufficient information to permit the design of a 
high-performance closed-loop system with the airplane as the controlled 
member. On the other hand, performance functions provide this necessary 
information in a form readily adapted to the problems of closed-loop design. 

This investigation is concerned only with some of the performance 
functions which describe the lateral motion of the airplane. 


Specifically, they are: 


(PF) al M,) 
(PF) 

a(é We] 
(PF) 

as. Wy) 


Pye .) 


The limitation in time available precluded the investigation of four 


more lateral motion performance functions of importance, ©.£. 


emia } 





PP C6 4B) 
(PYF ) ee | 
AST ) 


PY 
r 
Analytical expressions for these performance functions can be obtained from 
the differential equatione for lateral motion. 
The Laplace Transformations of the analytical equations are: 
(s - Yv + Ur - ef 2 0 


—vLT + (= 2 . Lr +(s- LP > “8 8 fe) 


i 
va. +(g- Nr + (- s- uP = A, b (se) + N. 6 (8) (1-1) 


“ae a r 
Side force and rolling moment due to rudder deflection and side forces 


due to rolling and yawing velocities are coneidered negligible. The axes 
and symbols conform to NACA notation. Fig. l-l pictures the axes. the 
symbols are defined in the Glossary, Appendix E. 

In order to measure the outputs due to a specific input, all other 
inputs must be zero. Letting > equal zero and specifying be as the 
forcing function, the solution of equation (1-1) yields the performance 


function, 


Pye W,) 
where Wo is equal tor 
This performance function for the B-253 ia computed in Appendix D, using 
values for the aerodynamic derivatives based on wind tunnel tests. Plots 
of the amplitude ratio and phase angle of this calculated performance 
function are given in Figs. 3 and 34, respectively. 


All of these lateral motion performance functions can be obtained 
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experimentally by sinusoidally oscillating the control surface and measur- 
ing the airplane’s steady etate response. This sinusoidal method requires 
a sine generator and an auto-pllot to supply the inputs, and requires as 
well, the usual measuring equipment. In addition, considerable flight time 
ie needed to obtain results over a sufficient range of forcing frequencies. 

In the case of rudder deflection inputs, a step-function input would 
yield a transient response from which the performance function could be 
determined. (See Ref .1) But a step-function of alleron deflection will 
produce a steady-state rolling velocity which in a short time will generate 
euch a large angle of bank that the motion could no longer be considered 
linear. 

The method of using a finite pulse of control deflection has the 
advantages associated with obtaining a transient response, and, since the 
control surface returns quickly to the trim position, the disadvantage of a 
etep-function input is eliminated. Finite pulses of elevator deflection 
have been successfully used to determine the airplane’s longitudinal 


(See Ref.2) But, to the knowledge of the authors, 


performance functions. 
this technique hase not been applied to the lateral performance functions of 
an airplane. 

The mathematical background for a method of obtaining the performance 
function from the transient response to a pulse is presented in Refs. (2) and 
(e). Briefly, the method consists of approximating the Fourier Transforms 
of the input and the response. The ratio of the latter to the former is the 


performance function, that 1s, 


(FT)[q,(t)) 
PP ta. ao] TMG, — 


or, more generally, 


at ia 





(LTLa., (t)] 


PP stay. ag) 7 GPa, Ce) |e = Sw (1-3) 


When the response ig convergent, it ia shown in (2) and (3) that 

—joAr w 

(FP)[q (t)]¥ (Fr)Cur(t)}., 2 a,(7)e s 
ea 


The lateral response of an airplane requires special treatment because 


(1-4) 


of the lightly-damped oscillatory mode and because of the divergent mode - 

found in many airplanes. The factored characteristic equation for the 

lateral motion of the B-25J, derived in Avpendix D, shows both characteristics: 
[ee + (2)(.241)(1.34)e + (1.34)"]Le + 2.9]Le - .0013)] = 0 (1-5) 
The lateral response of an airplane having this characteristic equation 

would not be convergent. However, by resolving the response into components, 

the performance function can be determined. 


For example, 


WPL (ont) Cape, )3 
- May tout? : (LT)(q,,‘* ] 
+ wr eifout foe 1) a (FPL 4 (out) (rem) §t? J (1-6 ) 
tan : J a tan t)] 


The Laplace Transform of Vout) (adv) st is obtained by first determining 





a= jw 


the analytical expression for Wout) (diy) **)s this is obtained by plotting 

(on semilog coordinates) the response, Vout) it) for large values of time, 

where, for practical purpoees, only the divergent mode remains. From this 

plot, the slope and intercept provide the necessary information to define 

the divergent mode analytically and the Laplace Transform of this analytical 
(ref. 8) 

expression can be determined. The Laplace Traneform of the input is 

equal to the Fourier Transform, since the input has the value of zero for all 

values of time less than zero and its integral is convergent. Thus, the first 


term on the right side of squation (1-6) can be found from a time response. 


No evidence of a divergent mode was found in this investigation for the 








B—-25J. Therefore, it was unnecessary to employ the technique described in 
the above paragraph. However, to illustrate the procedure the performance 
function of an idealized physical system with divergence was determined from 
a computed time response. This performance function was compared with the 
ideal performance function, thus showing the feasibility of the method. This 
analysis is contained in Appendix C. 

The divergent component of the time response must be subtracted from the 
recorded response. The integral of the remaining function of time converges. 
However, because of the lightly~damped oscillatory mode, the approximation 
of the Fourier Transform by triangles becomes tediously long. This problem 
can be overcome by subdtracting the oscillatory mode. The frequency, damping 
ratio, and amplitude of the oscillatory mode are determined as deacribed in 
detail in Apoendix D. This mode may be either subtracted over the region of 
time from zero to infinity, or it may be subtracted only in the region from 
any selected time, ty to infinity. In the first case, the performance 
function becomes: 


TG eat) Gace)? : C4 (out) (rem) 79) J FMC (UBT) (t)], 


PF & 
( L494] (PIgfuT(t)}, = q (7 eo 
mM 1.2.3 (in) ‘on 


~JjnAt w 
x (rr)(Ur(t)), M9 (cea ( tow) “Ta 
(FT)(UT(t)], Blac, .)(r,)Je es ae 
The unit triangle of t is here defined as a triangle of unit height and with 
a base length of c41T, as in Ref. 3. Similarly, the unit right triangle of 
t ie defined as a rignt triangle of unit height with a base length of Or. 
Note that T(out)(sew) is not the same function as T( eine in eq. 


(1-5), 


In the second of these two methods, the performance function becomes: 





FPN 2 cate) (ac)? lee $ ¢ S00) 


(PF) - ine 
( 14 (4n)4 (out) (FT)(OT(t)), zt a gaara 
Bat Die. « 


~ Jnr iw 
(FT)(UT(t)] | om ae ai nN 


n= Lager. 48 (17) 


(FT)(UT(t)], ze 
nee 192... 


If the Ort associated with (et ) equals that associated with Vin? only the 
summations remain in the second term on the right side of equation (1-7). 
However, if Or, is not equal to Or, a correction factor must de applied, 


since the Fourier Transform of a unit triangle of t is 


(sin we)” 
2 
(Fr) (UT(t)) = |——;—— ]_ 4 
(mer) 
2 

The Fourier Transforms of both a unit triangle and a unit right triangle, ag 
functions of wit, are plotted in Ref. (3). The Fourier Transforms of both a 
sine function and a cosine function over the region oe to infinity are 
developed in Appendix D, and their application to the determination of the 
lateral motion performance functions measured in thie investigation is 
discussed in detail. 

Either of the above methods is considered acceptable. However, the 
results of this investigation indicate that adding of the oscillatory mode 
at is ig a somewhat better method than including the oscillatory mode from 
time equal to rero. 

The input pulse should have sufficient area to produce a desirable 
magnitude of response, but should not cause motion beyond the linear region 
of the airplane's aerodynamic characteristics. Also, to simplify the data 
reduction, it is desirable that the pulse be of short duration. The pulses 


used in this investipation were roughly estimated from previous data on 


static sensitivities, such as p b/2U, Ref. 4, together with the experience 
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gained from the longitudinal pulee tests on this airplane. The linearity 
ranges of the linkages for the control deflection pick-offs (microsyn signal 
generators) and the sensitivities of the rate gyros were determined on the 
basis of those estimates with a 100 percent safety factor. After flight 
trials, the rough estimate of the aileron pulse was increased by about 

50 percent in magnitude, with the duration remaining as originally 
estimated. On the other hand, the estimated rudder pulse was reduced 


slightly after flight trials. 
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CHAPTER 2 


INSTRUMENTATION AND FLIGHT TEST PROCEDURE 


A. INSTRUMENTATION 


The airplane employed as a test vehicle for this investigation was a 


USAF North American B-25J, Army Number 44-30328. The airplane had been 


stripped of all armament. A photograph, a three-view drawing and a table 


of general specificatione and dimensions are presented in Appendix A. This 


airplane had previously been instrumented to record the normal acceleration 


and pitching velocity responses to a recorded pulse of elevator displacement. 


The amplifying and recording systems used in the earlier testing were retained 


for the present studies augmented by the following: 


dk, 


Adjustable rudder stop device (mounted on the co~pilot’s right 
rudder pedal), fulfilling the dual purpose of providing rudder 
pulses of any desired fixed magnitude, and of locking the rudder 
while aileron pulses were being applied. See Figs. 4 and 5, 
Appendix B. 


Adjustable aileron stop device (mounted on the co-pilot'’s control 
column), fulfilling the dual purpose of providing aileron pulses 

of any desired fixed magnitude, and of locking the ailerons while 
rudder pulses were being applied. See Figs. 1, 2 and 3, Appendix B. 


Microsyn signal generator units, linked to each aileron for recording 
the magnitude of the aileron displacement. See Figs. 6 and 7, 
Appendix B. 


Microsyn signal generator unit mounted at the rudder control dell 
crank (right rudder) for recording the magnitude of the rudder 
displacement. See Fig. 8, Appendix B. 


Rate gyro, for measuring the angular velocity in yaw. See Fig. 15, 
Appendix B. 


Rate gyro, for measuring the angular velocity in roll. See Fig. 15, 
Appendix B. 


Electronic current regulator, used to regulate a set value of etiffness 
current applied to the elastic restraint generators of the roll and yaw 


gyros. 


Four channels of amplification and four channels of a Consolidated 


recording oscillograph were utilized to record: 


a 


ee - an) a - 
is we oh 
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1. Aileron input pulses. 

2. Rudder input pulses. 

3. Roll angular velocity. 

4. Yaw angular velocity. 

The amplifiers and oscillograph are described in Appendix B. 

Several ground calibrations were made before any flights were 
attempted, and a final ground calibration was performed prior to the 
final flight. The methods of calibration and the calibration curves are 
presented in Appendix B. 

The center of gravity of the airplane was determined by means of a 
USAF electronic weighing kit. The center of gravity of the airplane was 
maintained near the center of gravity of the front main fuel tanks by 
shutting off the other tanks during flight. This caused the shift in the 
center of gravity during flight to be negligible. 

An accurate determination of the moments of inertia of the B-25J 
airplane for use in the theoretical calculations of this thesis posed a 
difficult problem. Obviously, they could be determined most accurately 
by experimental means, but, because of the size and weight of the B-2hJ, 
such a procedure was beyond the scope of the present study. As a practical 
alternative, the moments of inertia were fixed by careful estimates from 
known data already established on similar aircraft. 

Estimates were first based on the fairly complete data for a B-25J 
contained in Ref. 1. However, a subsequent comparison of the values for 
ey contained in the data with the value for ly determined for our test 
airplane by other investigators (Ref. 2) indicated that the test airplane 
and the reference airplane differed somewhat in configuration. Therefore, 
the final estimates were based on a comparison with similar data from the 


A-26 airplane discussed in Refs. 7 and 9. This comparison was expected to 
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increase the validity of the estimated figures because the test airplane 
and the A~26 had similar dimensions and differed only slightly in weight. 
Thies comparison yielded somewhat higher figures for “3 and 162 than those 
given in Ref. 1, in line with the higher walue for I yydetermined experiment- 
ally for the test airplane, to which reference has already been made. 
B. FLIGHT TEST METHODS 
Thirty minutes prior to a flight, all equipment was turned on 
temporarily and the following items were checked to insure proper 
operation: 
1. 300-volt @-c supply 
c. cé6-volt three-phase gyro wheel excitation 
3. 150-volt d-c supply 
4. Gyro stiffness current 
oO. Amplifier outputs 
6. Paper supply in oscillograph 
7. Frequency of the inverter output 
Soon after take-off, all equipment was turned on again = this time 
for the duration of the flight — to insure sufficient warm up time, and 
the items mentioned above were rechecked. Then, with the airplane trimmed 
at 175 mph and 10,000 ft. altitude, test pulses were applied so that the 
zero position of the light. traces of the oscillograph, the amplifier 
attenuation settings, and the input pulse amplitude could be adjusted. 
An effort was mde to keep the attenuation settings conatant so as to 
reduce the number of calibrations required. 
After the preliminary adjustments had been made, the airplane was 
carefully trimmed, with all controls held firmly against the stops. the 
pilot then started the oscillograph, and three seconds later gave the 


signal for the impulee. The oscillograph was stopped by the oscillograph 
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operator when the oscillations had attenuated to a negligible amplitude. 


The number of control stops used in these tests made it necessary for 


the pilot and co-pilot to follow a certain sequence of operations. This 


sequence was: 


a: 


Prior to recorded run} 


aA. 


Pilot trimmed airplane to maintain straight and level flight, 
holding the elevator control against the elevator stop. 


Compilot adjusted aileron and rudder stops to give desired pulse 
amplitudea. (This was done by applying pulses of various 
magnitudes, within the linear range, until the desired response 
was observed on the oscillograph. 


a recorded run: 


Pilot trimmed airplane in straight and level flight, holding the 
elevatox control against the elevator stop. 


Compilot moved adjustable trim screw for rudder until the proper 
stop was engaged, the rudder was then held firmly against this 
stop. 


Co—pilot moved the adjustable trim stop pin for the aileron until 
the proper stop was engaged, the aileron control was then held 
firmly against this stop. 


Pilot called oscillograph operator and, if the recording equip- 
ment was ready, pushed the button to start the oscillograph. 


Three seconds after starting the oscillograph, the pilot 
signalled the co-pilot to apply the pulse. 


On each flight, two records were made for each of the following types 


of pulse: 
ie 
e. 
3. 


4. 


Positive aileron (right aileron down) 


Negative aileron 


Right rudder 


Left rudder 


After the above records had been taken, the airplane was held ona 


straight couree at 10,000 ft. for a sufficient time to make calibration 


records. 


The method of calibration is discussed in detail in Appendix B. 
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Altitude 

Indicated airspeed 
Outside air temperature 
Fuel load 

Engine rpm 


Inverter frequency 


Gyro elastic restraint current 


Signal generator excitation current 
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CHAPTER 3 
RESULTS AND DISCUSSION 

A satisfactory instrumentation and calibration of the test vehicle 
was achieved. It is shown in Appendix B that the maximum uncertainty in 
calibration was 1.7 percent. 

The flight test data was reduced by the methods outlined in Ref. 3 
with the following results: ] 

(2) The uncertainty in reproducing results, obtained by analyzing 
three independent records of rolling velocity response to input aileron 
pulses, was founc to be 1.5 percent average deviation from the average 
measured performance function in amplitude ratio, and a average deviation 
in phase angle. (See Figures 3-1 and 32), 

(db) The comparison of the experimentally determined performance 
function with the performance function calculated from wind tunnel itability 
derivatives for rolling velocity response to aileron input showed very close 
agreement in phase angle at all frequencies within the range investigated 
(0.5 to 6.0 rad/sec). The comparison in AR showed both curves following 
the same general pattern with the calculated performance function having 
& lower AR at all frequencies except the resonant frequency. The resonant 
frequencies were displaced both in amplitude ratio and in frequency in 
accordance with the respective differences in damping ratio and undamped 
natural frequency of the oscillatory modes. (See Figures J~4 and 3-2). 

(c) The sensitivities of the calculated and experimental performance 
functions for rolling velocity response to aileron input compared as 
follows: 

5(A)(8_ W,)(cate) * ~*99? 2/eee 


Ss =~] .00 1 
(A) (8, ¥y) (exper) (aver) ‘ 





Utilizing vs te data obtained from flight test data (Ref. 4) an 


additional sensitivity is available, i.e., 


= -0.90 1/sec 
test data) 


(A) (6, HO Cre. 
(This sensitivity allows a rough check on overall calibration of the 
instrumentation eystem). 

(d) Comparison of the experimental performance function for yawing 
rate out to aileron input with the calculated performance function shows 
good correlation in AR. The PA comparison is satisfactory for low 
frequencies (.5 to 3 rada/sec). However, at high frequency a wide scattering 
of points is produced using three different systems of analysis. (See 
Figures 3-3 and S=4). 


(e) Sensitivity comparison of the calculated and experimental perform 


ance functions for yawing velocity response to aileron input is as follows: 


*(a)(6, MJeate = .c09 1/sec 


= .17 1/sec 
exper 


(a), Wo] 
(f) From an investigation of three different methods of dividing the 
output into comnonents, as discussed in Appendix D, it was found that the 
most accurate and farteat system of graphical analysis for use in the case 
of a-lightly damped oscillatory response was method (2) of Table 2, 
Appendix D. In this method the usual triangular approximation is used 
until the time response is only a pure damped sinusoid. At this point, 
C= 7 the Fourier Transform of a cosine or sine, whichever is appropriate, 
is added. 
(zg) Investigation of the scattering of points at high frequency in 
the performance function, FP) AD (8, > showed that the uncertainty was 


a function of the following: 
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(1) The measurement uncertainty in determining vector magnitudes from 
the oscillograph record. 

(2) Determination of the proper base Lane’. 

(3) The time response pattern — depends on whether the maximum ordinate 
occurs early or late in the time response, and on how large a proportion of 
the responee is purely oscillatory. 

(4) The ratio of magnitudes of the last plotted vector in the output 
summation to the resultant vector representing the response. 

(5) The method of component breakdown. 

(6) The frequency —— uncertainty increases with frequency in the usual 
case, 


(7) Size of At used in the approximation of the input and output. 


Instrumentation and the actual test flights for this project consumed a 
relatively large percentage of the time available. Consequently, the time re- 
maining for analyais of the oscillograph records required that only those 
analyses most pertinent to the objectives of the project be undertaken. Since 
the main objective of the project was to investigate the practical aspects of 
the application of the pulse technique to dynamic systems in general, and since 
determination of the performance functions of the test vehicle were of secondary 
importance, the scope of data analysis was limited to the following investigation: 

l. Determination of the uncertainty in repeating results by obtaining 
a performance function from each of three oscillograph records independently, 
using aileron pulse inputea and rolling velocity response. 

2c. Evaluation of several different graphical plotting methods as to 
accuracy, speed, and applicability to various response configurations. 

3. Comparison of the experimental performance function with the 


calculated performance functions, using response in rolling velocity and 


yawing velocity to aileron inputs as examples. 

4. Seeking out of limitations imposed on the accuracy of the pulse 
technique due to the approximations involved in the process of converting 
from the time domain to the frequency domain. 

Table 1 of Appendix D lists the oscillograph records analyzed, and the 
method of analysis used. | 

The uncertainty in reproducing results useing the pulee technique is 
considered to be entirely satisfactory, and well within the demands of usual 
engineering requirements. In evaluation of the ability of the pulse 
technique to reproduce results two investigations were pertinent. They 
were: (1) the ebility to produce the same performance function by analyzing 
different records, and (2) the ability of two different analysis methods to 
produce the same performance function when applied to the same oscillograph 
record. The results of the first investigation are shown in Figures &1 and 
S2. The three oscillograph records used to obtain the experimental results 
shown in Figures 8&1 and &2 are tabulated in Table 1, Appendix D. The 
experimental curve shown is the average of three performance functions deter- 
mined independently. The maximum AR deviation from the mean curve ig 5.9 


percent at w, = 1.5 rad/sec., and the average AR deviation from the mean 


f 
curve is 1.5 percent. For the phase angle, the maximm deviation is 6. 
and the average deviation is ie 

The uncertainty in repeating results by utilizing the cntes graphical 
methods listed in Table 2 of Appendix D, was found to be a variable depending 
on the response pattern in the time domain, and the forcing frequency. In 
general, the method described as Method (2) in Table 2, Appendix D, offers 
the least uncertainty. This method utilizes the same degree of approximation 


in determining the Fourier Transform of the input and output from the recorde 


in the time domain out to t = ° Disregarding freauency effects, the 
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uncertainty of this method may be held fairly constant with variation in 
response pattern by selecting % at increasingly higher walues as the 
response pattern becomes increasingly more oscillatory. Comparing the 
pattern of rolling velocity response with the pattern of yawing velocity 
response using aileron inputs for both, it was found that in the former 

a low proportion of the response was oscilletory while in the latter a high 
proportion of the response was oscillatory. Consequently, in the latter 
case, if nn was not increased beyond the ‘ used in the rolling response 
case, a higher proportion of the response was represented by a Fourier 
Transform which is exact, while the input remains an approximation to the 
exact Fourier Transform. The ultimate limit is to let o increase to a 
time beyond which the response is negligible, plotting all vectore directly 
fron the record. This is listed in Table 2, Appendix D as method (1). 

Thiea method offers the least error in the hypothetical case where 
instrumentation inaccuracies and noise are at such a low level that 
ordinates may be measured accurately at large values of time. 

Increasing ” as the percentage of oscillatory portion of the response 
increases is also important in nullifying errors occassioned by choosing . 
at a point which is not exactly a true peak or zero of the response curve. 
An error of this nature serves to shift the phase angle of the added vector 
for the oscillatory remainder (t 2 ,) vy an amount equal to w At. If the 
magnitude of the last vector plotted in the summation of triangles is only 
& small percentage of the added vector, the error in phase angle may be 
large. However, if the magnitude of the vector added is of exactly the eame 
length as the last vector of the summation, an ornor in ta as large as t = Ar 
will yield no error in the PA of the response vector. 

The method of dividing the output into components, described in Appendix 


D, Table 2, aa method (3), appears to be subdject to the greatest inaccuracies 


co" 








any of the methods investigated in analyzing a lightly damped oscillatory 
response. Theoretically, this method of analyeis is as sound as any of the 
methods investigated. Howewer, the practical application of thir method 
presents an additional source of error since the input and output are not 
represented with the same degree of approximation. In the lightly damped, 
highly oscillatory response, essentially the entire output is represented 

by the Fourier Transform of the oscillatory mode, which is the exact 
transform of an analytical expreesion. On the other hand, the input is 
approximated by a vector addition, the exactness of which depends on the 
accuracy of curve fitting achieved. At relatively high frequency, the 

emill resultant vector of the remainder summation, although approximated 

with the same degree of exactnees as the inout resultant vector, is 
ineffectual in changing the total output phase angle which is dominated by 
the large vector representing the oscillatory mode. Consequently, the output 
phase angle becomes nearly static as the limiting phase angle of the oscillatory 
mode is approached, while the input phase angle continues to grow as the 
forcing frequency increases. 

In performing the component breakdown necessary in this system, additional 
measurement errors may accumulate if more than one set of measurements are 
taken from the response. That is, in removing the oscillatory mode, the 
ordinates ae determined by the analytic expression of the oscillatory mode 
should be subtracted algebraically rather than graphically from measured 
ordinates of the total response. Furthermore, since the oscillatory mode 
represents the major part of the response, great care must be exercised in 
determining Wee DR, and esi: 

In Figures 33 and 3-4, the results of analyses by the various methods 
discussed are shown. The scatter of points at high frequency shows the 
inaccuracies of the several methods. Theoretically, the PA curve should 


level off in the vicinity of -90° phase angle. However, Fig. 3-4 shows the 
—2o" 





experimentally determined PA lagging less at w, of 5 and 6 rad/sec than 


r 
at 3 rad/sec. In investigating this discrepancy between theoretical and 
experimental results, the question of accuracy in approximating the Fourier 
Transform arose. 

To determine how closely the approximation of input and output time 
response by summation of triangles fits the actual curves, frequency spectra 
of yaw response and aileron input were prepared by numerical means. The 
results are shown in Figures 35 and &6. The circled points show the 
results of the numerical analysis and the iotibetl Gaeeb- Gilets the mean 
curve of the Aetbnenied points obtained in Figures 33 and 24. The phase 
angles determined by the numerical analysis show a variation of more than 30° 
from the average curve at many points. The explanation for such a discrepancy 
wae found to lie in a difference in measurement of the ordinates, (q(t), 
between those used in the numerical analysis and those used in the graphical 


analysie. The measurements of the respective areas, for each analysis, are 


as follow: 


AREA AREA 
Numerical Graphical 
Analysis Analysis 
i positive lobe 0.268 in sec O.c6 in sec 
t 
1°” negative lobe 0.6170 in #ec 0.5990 in sec 
nd 
2 positive lobe 0.1126 in sec 0.1096 in sec 


In the two most important lobes, the deviation in measurement of srea 
is approximately 3 percent, and in each case the graphical analysis is 
the smaller area. Thus it is noted in Figure 37 that a emall uncertainty 


of 3 percent in area determination may cause as much as 30° difference in 


4 


phase angle, provided, the vector plot of the output ‘'winds up’’ near the 


origin. Since the numerical analysis gives only a check on how closely 


the triangular summation approximates the area of the curve, the performance 
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function was again determined at several forcing frequencies using the same 
ordinates as used in the numerical analysis. The triangles in Figure 3-6 

show the points thus determined. It 16 seen that the uncertainty of repro- 
ducing results is satisfactory, indicating that the use of triangles in 
approximating the area yields approximately the same accuracy as the numerical 
approximation. It is now apparent that any inaccuracies in measuring the 
ordinates hweacritical effect on the performance function when the end 

point of the vector summation is near the origin. However, the PA variations 
due to inaccuracies in measurement of ordinates is not constant, but is a 
function of frequency. Although the measuring error in this case gives an 
increase of ordinate magnitude in both the positive and negative lobes, the 
errors in a positive lobe do not nullify the errors ina negative lobe 

except at certain frequencies which are indeterminate unless the actual 
measurement errors are known. Figure 37 shows the resulte of measurement 
error at vw, = 5 rad/sec. Due to the phase angle shift with each added 

vector (w Ar) the errors are apportioned in accordance with the sine and 
cosine of the vector phase angle multiplied by the magnitude of the error 

in that vector. If the sum of the imaginaries and reals of the measurement 
errors add up to give a vector which has an angle corresponding to that of the 
original resultant vector, no change in PA will result. However, this rarely 
happens, and at high frequencies a large phase angle change results. Note in 
Fig. &7 that the largest differences in the vectors of the first two lobes 
occur in measuring vectors which have phase angles close to +90, i.e., vector 
10 in the first lobe and vectore 28 and 30 in the second love. All of these 
difference tend to increase the phase angle at We = 6. Bach frequency will 
have a different total PA error due to measurement uncertainties, since the 


phase angle of each vector changes with each different forcing frequency. 





Again note that measurement errors have an insignificant effect unless 
the et oni of the summation is near the origin. At a, 1.75 rad/sec 
the PA is the game, whether determined graphically or numerically, and the 
effect of measurement differences on uncertainty in determining the 
performance function is negligible as shown by the superimposed points 
lying on the average experimental curve. 

In general, it appears that any uncertainty, whether in measurement 
of ordinates or base line location, gives a magnified uncertainty in PA 
if the end point of the summtion ia near the origin. 

The phase angle uncertainty is not due to a breakiown in the theory of 
the pulse technique, nor caused directly by the triangular approximation 
of the input and reeponse. Rather it ie a function of uncertainty in 
measurement of ordinates and in plotting which accum late under the 
conditions of high frequency (rapid phase angle change of each added 
vector) and low amplitude ratio (requires and point to be near the origin). 
Furthermore, this uncertainty in phase angle is not restricted to the 
predominately oscillatory response. It merely appears at a lower frequency 
in this case due to the rapid attenuation in amplitude ratio (small 
vectors in the early portion of the response) . The same phase angle 
uncertainty would occur in the rolling velocity reponse to aileron 
input if the performance function determination had been carried out to 


those high frequencies where the amplitude ratio was greatly attenuated. 
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CHAPTER 4 
CONCLUSIONS AND RECOMMENDATIONS 
A. CONCLUSIONS 

Ase a result of this investigation the following conclusions are 
submitted: 

1. The pulse technique is a practical and suitably accurate method 
for determining aircraft performance functions from flight tests. 

@. Economy in flight time and inetrumentation are gained by the 
pulse technique as compared to the einusoidal method. 

3. The accuracy of the resulte obtained is enhanced by the short 
duration of test flights in that the changes in mass and moments of 
inertia due to fuel consumption are negligible. 

4. Erroneous performance function phase angles will in general be 
obtained at relatively high frequencies when the amplitude ratio has 
attenuated to a very small value (order of five percent). This is due to 
the fact that the errors in measurements of the ordinates of the recorded 
reeponse accumulate at those frequencies where the series of vectors, 
representing the ordinates and their individual phase angles, encircle 
wne*oripin. 

5. Lightly damped oscillatory responses may be analyzed with reasonable 
simplicity and accuracy. However, in the case of a response which is almost 
wholly oscillatory, uncertainty in the data may cause anomalous phase effects 
to appear in the plots of the performance function at the higher forcing 
frequencies when the amplitude ratio is considerably attenuated. These phase 
effects are more pronounced in the performance functions which are highly 
dependent on coupling terms in the equations of motion rather than the 
forcing function itself, since high attenuation in AR occurs at a relatively 


low frequency. For example, FF and PF r 
q 1 ce Va) [ } are such performance 
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functions, 
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6. <A predominately oscillatory response is best approximated by use 
of isosceles triangles to a value of time beyond which only a dateped sinusoid 
remains. The Fourier Transform of the sinusoid is then obtained by trens- 
forming the analytic expression which is then added vectorially to the 
approximete transform of the early portion of the response. 

7. The most important sources of error in reduction of data are: 

a. Choosing the base line on the oscillograph record. 
bd. Measurement of the ordinates from this base line. 

8. The calculated and experimental performance functions are 
qualitatively comparable, but the quantitative differences indicate the 
requirement that the dynamic charecteristics of all high speed and non- 
conventional aircraft be determined by flight tests. 

B. KECOMMENDATIONS 

To aid further research in the use of the pulse technique to determine 
aircraft performance functions the following recommendations are presented: 

1. Obtain relatively long records of lightly damped oscillatory 
responses to facilitate the determination of We and DR of the oscillation. 

2. Amplify low amplitude predominately oscillatory responses as mich 
as is consistent with the noise leval and linearity of the recording 
instruments, the friction torque level of the gyros, and the linearity of 
aircraft motion. 

3. Limit the shock mounting of the gyros to landing and take-off 
only. During tests the gyros should be rigidly attached to the aircraft 
structure so that the input to the gyros is the true angular velocity of 
the airplane and is not modified by the flexibility of the shock mounts. 

4. Superimpose a reference trace on each response to aid in the 


determination of the base line and in indicating any drift or divergence 





of the response trace. 

5. Investigate by flight teet the lateral motion performance 
functions of an airplane having a relatively high degree of epiral 
divergence. Since the degree of spiral divergence increases with angle 
of attack, it may be necessary to perform these measurements at a high 
angle of attack in order to obtain suitable responses. 

6. Investigate the performance functions of an airplane having 
dynamics quite different from those of the B-25J. Several nodern fighter 
typee exhibit an undamped oecillation about the Z axie. A study of this 


motion ehould prove most worthwhile and interesting. 





APPENDIX A 


DESCRIPTION OF AIRPLANE IN TEST CONFIGURATION 


The airplane used in the tests was a USAF B-25J under bailment contract 


to the Instrumentation Laboratory, M.I.T. The airplane had previously been 


stripped of all guns and turrets, as shown in Fig. A-l. 


Following is a listing of the specifications and dimensions of the air- 


plane as it was used in the tests described in this thesis: 


Manufacturer 


Type 
Serial Number 
Center-of-gravity range 
Forward limit 
Aft limit 
For tests 
Gear retraction moment/1000 
Overall length 
Height 
Wing, airfoil section 
Root 
Tip 
Span 
Area (Total) 
Area (less ailerons) 
Aspect ratio 
Taper ratio (tip chord/root chord) 
Chord at root 


Chord near tip (380.75 from fuselage) 
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North American Aviation 
Corporation, Inc. 


B=2OU—25- CC . 


44-30328 


16.5 percent M.A.C. 
34.0 percent M.A.C. 
27.4 percent M.A.C. 
59.7 in-lbds 

Bo ft 8.75 in 


16 ft 4.19 in 


NACA 23017 
NACA 4409 R 

67 ft 6.704 in 
609.8 aq ft 
577 167 #q T% 
7.48 

0.415 

le ¥% 1076 in 


6 ft 10.14 in 





Projected tip chord 
Mean aerodynamic chord 
Length 


Distance of leading edge back of 
nose reference datum line 


Incidence 
Root 
Tip 

Dihedral (25 percent line) 
Wing center section 
Wing outer section 


Sweepback (leading edge) 


Aileron 


Wing 


Type 
Area 
Span 
Chord 
Inboard 
Outboard 
Travel 


Up (from neutral) 
Down (from neutral) 


Tab area 
Tab span 
Tab chord (mean) 


Tab travel 
Up (from aileron trailing edge) 


Down (from aileron trailing edge) 


flaps 


Type 


Area (total) 
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o ft 4.257 in 


116.16 in 


212.68 in 


SO, OF 


QO 2" Bor” 


4 28° pa’ ' 
QO “12° %8"" 


- let ey 


Internal balance unsealed 
32.13 aq ft 
137 .25 in 


19.87 in 
12.45 in 


15° 
3.86 aq ft 
57.911 in 
5.070 10 
12° 


12 


Variable-slot 
trailing edge 


75.8 sq ft 





Center section flape 
Area 
Span (mean) 
Travel (down) 
Outboard flap 
Area 
Span (mean) 
Travel (down) 
Horizontal tail surface 
Area (including elevators) 
Span 
Maximum chord 
Incidence 
Dihedral 
Elevator area (two, including tabs) 
Elevator span (two) 
Elevator travel 
Up (from streamline 
with stabilizer) 
Down (from streamline 
with stabilizer) 
Elevator trim tab area (total) 
Elevator trim tab span (2) 
Elevator trim tab mean chord 
Elevator trim tab travel 
Up (from elevator trailing edge) 
Down (from elevator trailing edge) 
Vertical tail surfaces 
Area (total) 


Yin area (two) 


Span (maximum) 


32.6 oq ft 
127.5 in 


43° 


43.2 eq ft 
cue. in 


46° 


132.4 aq ft 
ee ft 2 ia 

7 ft 1.875 in 
2° 

None 


50.6 sq ft 


216.5 in 


26° 


10° 


4.22 eq ft 
S7.2 


7.875 in 
20° 
91 eq ft 


47.8 eq ft 


104.25 in 





Rudder area (two including tabs) 43.c aq ft 
Rudder span (maximum) 101.1 in 
Rudder travel 


Right (from etreamline with fin) 20 
Left (from streamline with fin) 20 


0 
0 


Rudder trim tab area (total) 3.18 sq ft 
Rudder trim tab span 36.4 in 
Rudder trim tab mean chord 6.3 in 


Rudder trim tad travel - 
Right (from rudder trailing edge) le) 
Left (from rudder trailing edge) 12 


Fuselage 
Maximum width 4 ft 8.5 in 
Maximum height 7 ft 4 in 


Length (tip of nose to tip of tail) 63 ft 5.75 in 


Engine 
Humber 2 
Type Wright Cyclone 
Designation Re2600-13 or —29 
i 
Fumber of cylinders 14 - double row 
Gear ratio 16.9 
Supercharger gear ratio 
ow 7 06:1 
High 10 .06:1 
Power settings: 

Condition rpm Manifold pressure (In of H ) 
Take-off 2600 46 . 
NEP 2400 40 
Cruise 2100 29 
Approach 2000 20 
For teste 1800 24 
Propellers 

Manufacturer Hamilton Standard 
Type Controllable pitch 


~40- 





Rumber of blades 
Diameter 
Blade type 


Control type 


Hub type 


Pitch et 
Low (fine 
High (coarse) 


Landing gear (estimated data) 


Main gear assembly 
Weight of retractable portion 
Radius of gyration of retractable 
portion about retraction pivot 


Nose gear assembly 
Weight of retractable portion 


Radius of gyration of retractable 
portion about retraction pivot 


Weight 
Moments of inertia (estimated) 


yx 


ty 


yz 


120 ,000 


3 
le ft 7 in 
6359 A-18 


Hydromatic full 
feathering 


2SE50-473 


0 


a 
90 


600 1b (each) 


48 in 


200 lb 
46 in 
26,637 1b 
63,000 slug rt” 
slug rt” 


2 


-1930 slug ft 
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APPENDIX B 
INSTRUMENTATION AND CALIBRATION 

A. INSTRUMENTATION 

As noted in the Introduction, finding the performance function by means 
of a pulse input has the advantages of keeping the response within the 
linearity limits of deviation from the equilibrium flight conditions, and of 
requiring a minimum of inetrumentation. The devices required for applying 
the input pulses were designed to accomplish the following functions: 


1. Allow rotation of the control surfaces between any desired limits, 
with positive stops at the ends of the arc of travel. 


2. Allow variations in the size of the pulse. 

3. Allow the pulse to be applied in both directions, that is, allow 
changes in the position of the mechanism relative to the pilot's 
controls without altering the pulse size. This function also 
allows the pilot to hold two controls rigidly against stops while 
the aircraft is responding to a pulse imparted by the third control. 

4. Allow quick disconnect so as to make it possible to remove all 
impediments from full use of the pilot's controls in case of 
emergency. 

The mechanisms employed were chosen from many possible designs ae the 
simplest that could perform the above functions and still be fabricated with 
a minimm of machining. The aileron puleing mechanism, shown in Figures 5-1, 
B-2, and B-3, coneiste of a plate mounted on the right control wheel which is 
slotted to receive two pivoted bars which are set to give the desired size of 
pulse. <A slotted length of duralumin angle, hinged at one edge, was mounted 
on a shelf fastened to the top of the control column. The elot allowed a 
single stop to be positioned so that the chosen pulse could be initiated at 
any trim condition by moving the single stop against either one of the stops 
on the wheel plate. The hinged edge allows the etop to be thrown out of 
position in order to clear the aileron control in case of emergency in flight. 


Approximately 45 degrees of wheel rotation were used in applying the desired 


aileron pulse. 





The mechanism used for applying rudder pulses is shown in Figures 
B-4 and B-5. Choice of the design selected was dictated by the limited 
installation space and by the requirement that the pulse setting and quick 
disconnect units mist be easily accessible to the co-pilot in flight. The 
mechaniem in its final configuration consists of a long rod, one end of which 
is attached to the right rudder pedal. The other end of the rod slides ina 
block which isa pivoted at the center, and mounted on a plate fixed to the 
aircraft structure. The block receiving the sliding rod is cut away at the 
top so that a stop pin welded to the sliding rod may move in the slot the 
provided. The cover of the block is made up of a frame hinged at the inboard 
edge and containing two stop bars connected by a positioning screw to provide 
adjustment of pulse size. In addition, both bars, holding the set pulse size, 
may be moved by means of a second screw to permit engaging the pin on the 
sliding rod with either of the stop bars when in the trim condition. The 
hinge on the inboard side allows the block cover to be moved clear of the rod 
stop, thus quickly and positively freeing the rudder controls. Figures B=-4 
and B-5 show the cover block clear of the rod stop in the stowed condition. 
Approximately G43 in. of linear throw of the rudder pedals provided the 
desired rudder pulse size. 

The microsyn pickoff units inetalled for measurement of control surface 
rotation were of the variable transformer type designed by Dr. B.K. Mueller 
of the M.I.T. Instrumentation Laboratory. The installation of these units 
in the right and left wing are shown in Figures B-6 and B-7 respectively. 

The mounting positions and linkage lengths are identical in each wing. The 
final determination of linkage seneitivity, that is, piokeff rotation per 
control surface rotation, was governed by accessibility in the wing or rudder 
cavity, linearity limits of the pickoff (approximately 150 milliradians) ,and the 


maximum control surface motion anticipatec in obdtaining the desired magnitude 
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of response. Shielded pickoff leads were used throughout, and, as shown 

in Figs. B-6 and B-7, safety bracketea were attached to the microsyn mountings 
to prevent any possible damage or jamming of controle in case the linkage 
worked free at the pickoff rotor. 

The rudder pickoff installation is shown in Fig. B-8. The linkage is 
similar to those ured in the wings, the bell crank at the rudder control 
serving as the initial point of motion measurement. 

The rate gyros employed were elevation and deflection gyros developed 
by the Instrumentation Laboratory, M.I.T. The characteristice of these 
gyros were modified to suit the measurement requirements of the project by 
changing the damping, and were installed as roll and yaw rate gyros 
respectively. 

Laboratory calibration of the gyros prior to installation in the aircraft 
was necessary to determine the excitation currents desired for the elastic 
reetraints and the signal generatore (pick-offs), the temperature required 
for the damping fluid, and the friction torque level. 

The friction torque level was first determined to judge whether the gyros 
would provide a satiefactorily sensitive response to allow measurements of 
the desired accuracy and precision. The final results are shown in Figures 9 
and 10 for the roll and yaw rate gyros, respectively. The maximum friction 
torque level for the roll gyro was 50 dy cm, and for the yaw gyro was 140 
dy cms, and the unbalance was reduced to a torque level below that of the 
friction level. The maximum friction torque level for thé roll rate gyro 
represents an angular velocity input of 0.25 milli-radians/seec and for the yaw 
rate gyro, the maximum friction amet represents an angular inovut velocity of 
O.7 milli-radiana/sec. This friction level is connidered satisfactory since 
the oscillograph record ordinates can be measured only to a maximum accuracy of 


Ol inch, the latter figure representing approximately 1/2 milli-radians/sec in 
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most of the recorde. 
The rate gyro equations of motion were investigated to make certain 
that the desired characteristic could be obtained. The general equation 


representing both gyros is: 


A + C. A + 8 


"C) *8(rtnx) 4 M,) 7 "(e) (eet) (er) (a.m_) 4 


where the subscript c refers the indicated quantity to the computer shaft, 
and the subscript g refers the indicated quantity to the gimbal shaft. 
This equation may be reduced to the general form of the second-order 


differential equation yielding 





A 
ao 5. ne, 
Hi) €e) (WA) Oy (yey?) + 5 20k) (FR) 
where 
(DR) “8 
DR) = 
ot (arr) mr 
ey tee, ya S(er) (4M) _ 
a 1 (eff) 


fo determine the gyro undamped natural frequency, the effective moment 
of inertia at the computer shaft was measured using low elastic restraint 
currents and no damping. The period was recorded on a brush type oscillograph, 
and the moment of inertia was then determined from the following equation: 


Ss 
(ex) (aM) r = 264 ‘Fumite- ber 


Terr) * aa 
Rechecks at different values of elastic restraint verified the above 

figure. To establieh the amount of elastic restraint desired on the gyros, 

maximim rates of roll and yaw for the magnitude of control eurface deflection 


inputs anticipated were determined from B-25 roll rate data (Ee) and fron 


Cornell Laboratory Reports covering fate responses to step inputs of control 
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deflection (Refs. 4 and 5). The maximum rates determined were 0.4 rad/sec 
in roll and 0.5 rad/sec in yaw. The elastic reetraint unite were then 
calibrated for both gyros, as shown in Figs. Bell and B-12, using various 
values of elastic restraint current. The maximum rate tabulated for each 
sensitivity indicates the rate input for the particular sensitivity which 
will rotate the gyro computer shaft 150 milliradiane. At deflections 
exceeding this value, the signal generators become nonlinear. The overall) 
sensitivities of the elastic restraints, considering both excitation current 


and angular displacement as inpute are: 


2 = 1820 Syne ch 


; mM) (ma)? rad 


"reantwe) (ay @ 


* 612 Syne oR 
(ma)~ rad 


* (yaw) (er) (A, é M) 
Selection of the desired undamped natural frequency depends on the 
desired frequency response of the gyro. The rate gyros mst accurately 
reproduce aircraft response below about 10 rads/sec, offering negligible 
or minimum attenuation and phase shift, while at the resonant structural 
vibration frequencies of the airplane, the attenuation mest be a maximun. 
Previous instrumentation work on the B-25 has indicated that the most 
pronounced resonant frequencies were found at 16 and 30 cps. To fulfill 
the response requirements outlined above, the following acceptable limits 
for the gyro characteristics were determired: 
Undamped Natural Frequency, — & - 8 eps 
Damping Ratio, (DE) = 0.6 - 1.0 
With the above ranges, elastic restraint currents were selected to satisfy 
both undamped natural frequency and maximum rate input requirements. 


A damping fluid was then selected to fit the above requirements of DR. 
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Figure B-13 shows the viscosity vs temperature characteristics of damping 
fluid DCC 350, which was the fluid used in both gyros. DCC 250 offered the 
amallest gradient of viscosity to temperature of any of the fluids available 
having the viscosity range desired. 

The following chracteristice, coefficients, and sensitivities describe 


the gyroe as used in the final teat flights: 


ate Yaw rate gyro 
2 2 
1.) (err) 264 dy cm sec 264 dy cm sec 
dy cm cx 
“(c) 17200 ua rad/ sec 
Damper temperature 160° r . 177° ¥ 
So dy cm 5 dy cm 
8 (or) (AM) 3.6 X 10 cod 4.665 X 10 oa 
v, 38 rad/sec 42 rad/sec 
a : 6.05 cps 6.7 cps 
(DE) 857 1.08 


The gyros were mounted in a dustproof container and secured in the air 
plane on an adjustable mount which allowed aligning the gyros with the 
horizontal. Since the B-25 trims at 175 mph in a nose high attitude with 
its horizontal reference line 3.6° above the horizontal, the gyro container 
mitt be tilted accordingly. The adjustable mount is shown in Fig. B14, 
and the gyro conthinen in position on the mount is shown in Fig. B-15. The 
gyro container utilized three thermostatically controlled heaters to maintain 
a fixed ambient temperature for the gyros. 

Figure B-16 shows the amplifier panel. Four channels of amplification 
were used to amplify the pick-off voltages representing total aileron deflection. 
rudder deflection, angular velocity in yaw, and angular velocity in roll. 

Each amplifier channel performed three functions (1) amplified the AC pick-off 
signal, (2) demodulated the AC signal, and (3) amplified the resulting DC 


signal. The output of each amplifier channel is controlled by 12 attenuation 
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settings in order to adjust oscillograph current to give the desired 
oscillograph displacement. 

The oscillograph used is shown in Fig. B-17. It is a 12-channel 
Consolidated Recording Oscillograph. The dynamics of the oscillograph 
galvanometer units have a negligible distortion effect in reproducing the 
inputs introduced by the amplifier outputs. The natural frequency of the 
units used is 120 cps and the DE is0.7 approx. 

The stiffness current regulator shown mounted in position in Fig. Bel6 
is a closed-loop feedback regulator functioning to hold the stiffness current 
at a set current within 1 percent wariation. Figures B-15 and B-17 show the 
exposed regulator chassis, the chassis having been removed from its case. 

To integrate the added instrumentation with the instrumentation 
previously existing in the aircraft, an additional junction box was 
installed. It’s location is noted in Fig. B-17. 

The three-phase, 400-cycle thyratron controlled inverter initially 
installed was found to give unsatiefactory regulation for the needs of 
this project. Consequently, it was replaced by a USAF PU-16 single-phase 
inverter, which was found to give satisfactory regulation. This change 
necessitated use of a phase splitter and stepdown transformer to obtain 
the 3 phase, c8v gyro wheel supply. A block diagram of the overall 
inetrumentation system is shown in Figure B26. 

As shown in Fig. B-c6 the output of each pickoff was loaded with a 
1000—ohm resistance to reduce the noise level and the null voltage. A 
variable resistance and a variable inductance were placed in series with 
the pickoff excitation coils (P.O. coile in series). The resistance was 
used to set the desired excitation current, and the inductance was used to 


bring the pickoff outputs into proper phase. 








B. CALIBRATION 


The objective of the measuring system calibration for this project 
was to establish overall sensitivities, expressing inputs to the measuring 
instruments in terms of output linear displacements at the oscillographe. 
The sensitivities are expressed in such a manner that multiplication of any 
displacement on the oecillograph record by the appropriate sensitivity yields 


the measured quantity in the desired units. These sensitivities are: 


S degfin 
(ail) (xX. . é) 

S(rua) (x 8) deg/in 

osc or 

5 depg/sec 
(roll gyro)(X W) in 

5 deg/sec 
(yaw gyro) (x Ww) in 


For purposes of simplicity and system flexibility, the overall sensitivities 
listed above were broken into two component sensitivities, one representing 
the instrument sensitivity and one representing the indicating system 
sensitivity. The inetrument sensitivity is a fixed quantity determined in 
the case of the control surface rotations by the linkage sensitivities and 
the signal generator excitation current, and in the case of the rate gyros 
by the elastic restraint current and the signal generator excitation current. 

The indicating system sensitivity may be varied at will by changing the 
attenuation setting at the amplifier. The amplifier sensitivity allows the 
desired size of recorded input pulses and output anguler velocities to be 
selected. Selection of the proper size of recorded displecements is governed 
on one hand by the linearity limits of the oscillograph trace and on the 
other hand by the eize record required for ease of analyeis. 


In order to obviate the necessity for calibration of each amplifier at 
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each attenuation setting, a ‘’master’’ pickoff and a calibration voltage 
source, the latter providing l2 voltage stops, were employed. The master 


pickoff was calibrated (Figure B-18) to determine the limits of linearity. 


m volts 


divieion ' “'™ 3 


The sensitivity determined was 5 » = 2 


mpo) (A rae 


master ~p.o. divisione as the limit of the linearity range. The mater 
pickoff was then used to "'buck-out’’ or neutralize the voltage produced 

by the instrument pickeffs. Thue, inetrument sensitivities were obtained 

doth in terms of voltage output for input rotation or angular velocity and 

in terms of _ pickoff output angle for input rotation or angular velocity. 

The master piokeff is shown in Fig. B-19. The master pickoff could be 
plugged into any of the pickoff outputs by means of bayonet jack recepticles 
on the amplifier panel. As shown in Fig. B-26, the particular pickoff out- 
put could then be switched from the amplifiers to the master pickoff. 

The procedure for obtaining the instrument sensitivities was straight- 
forward. For the ailerons, a propeller bubble protractor, measuring to 0.2 
degree accuracy was used to measure the angle of both ailerons with respect 
to the horizontal, and the total aileron differential angle thus obtained. 
For each measured oe the corresponding voltage was recorded. This voltage 
was then nulled by the MPO and the corresponding angle of the MPO recorded. 
The voltage representing total aileron angle was achieved by adding the 
voltages of the two aileron pickoff outputs in series. 

For calibration of the rudder deflection to MPO angle, the same system 
was used as outlined for the ailerons with the exception of the means of 
measuring rudder deflection. Figure B-20 shows the method employed. The 

center of rotation of the rudder was determined by means of a plumb bob. 
A piece of mirror was then attached to the rudder at the center of rotation, 
| and the plumb bob lined up with the center of a sector of a circle, the 


| latter being graduated in mile. Measurements were then read by means of a 
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: telescope equipped with a hair line. 





Figure B-2)1 shows the set-up used to calibrate the gyros. On a small 





turntable, having a selection of 4 speeds, an inclined plate was mounted to 
receive the gyro container. The can was 80 aligned that the spin axis of 
the roll rate gyro was rotated a measured angle about the gimbal aris. With 
the two gyros properly aligned within the can, this also rotated the input 
axis of the yaw gyro about its spin axis by the same measured angle. Asa 


consequence, the input to each gyro was calculated in the following manner: 


Roll gyro Wy a Wan sin @ 
Yaw gyro x, = Ving cos 8 


where War le the angular velocity 
of the turntable, and @ is the meaeured angle of tilt. 

The calibration curves for all four measuring instruments showing the 
measured quantity vs MPO divisions (each div. = approx. 2 mils) and, for 
the gyros, the measured quantity vs pickoff output voltage are shown in 
71e@. S-ce through 8-25. 

The instrument sensitivities thus determined are 


‘ deg 
(aileron] — div 


S(at1)(e, 8.) °(MPO)(Aypo €5) ~ 7(ai1)(M) (Ayo 8) 





: ; deg 
S(rud)(e,. 8.) S(MPO) (Ay 0.) ~ Crud) (M) yay 6.) ~ °°P8° aay 


= 0.994 mils/sec 
i di 


vy MPO 


( rudder] 


[roll gyro} (re) (e, Wy) (MPO) (Ayo e.) 7 (rg) (M) (Aypo Wy) 


ils/sec 
[yaw gyro] S S = § = 1.228 mils/sec 
Vv 
(ye)(e, W2) (MPO) (Ap 85) (ye) (¥) (Aung Wy) div MPO 
The indicating system sensitivities were determined on each flight 
whenever data was recorded. The indicating system sensitivity is made up 
of the amplifier sensitivity, oscillograph sensitivity and master pickoff 


sensitivity, that is: 





*(ind eye) (X46: Aypo? = S(upo) (6 Aypo) * (amp 1 ,¢) (ose) (X54. ‘) 


Upon the completion of each series of flight rune the measuring 
instruments were switched away from the inputs to the amplifiers, and a 
l1&step voltage source ewitched onto the amplifier inputs (shown in Fig .B-19). 
As the voltage was varied through the 12 eateps, a record was made on the 
oscillograph. These same le voltage steps were then applied to the master 


pickoff, the latter instrument being rotated until each calibration voltage 





was nilled and the corresponding output angle of the master pickoff was 
recorded. Thus increments of oscillograph displacement corresponding to in- 
crements of MPO angle were available for each teat flight, and an overall 
indicating system sensitivity was available without knowledge of the 
oscillograph and amplifier component sensitivities. That is, 
o Aveo 
5 (ind eys)(X. A oo) ~ @ x ec 
The overall system sensitivity for each measuring channel was thus 
available by combining instrument and indicating system sensitivities. 


Therefore, the sl expression of the sensitivity is 
8(q chan) (X00 a) ~ §( )(M)(dypoe a) °C (4nd sys) (X56. Ayao) 


To determine the system errors and uncertainties, etatic teste were 
made on the ground. A kmown value of each quantity to be measured was 
applied to the respective measuring instruments, and an oscillograph 
record made of the outputs. The calibration source voltages were then 
recorded at the same attenuation settings, and a record of master pickoff 
angles recorded for each calibration voltage. In this manner the overall 
sensitivity could be compared with the product of the instrument sensitivity 


and the indicating system senaitivity. 
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The linearity limite and of the measuring system components are as 


follows: 
Oomponen: Linearity Range 
Master P.O, +300 dive (150 mils) 
Ajileron Defl. 30° total deflection 
Rudder Def). +7° rudder deflection 
Roll gyro +0.3 rad/sec 
Yaw gyro +0.35 rad/sec 


Oscillograph defl. +2¢.7 in (approx.) 


The largest source of error in the measuring system lies in the 
insensitivity of the voltmeter used as an indicator when nulling a voltage 
with the master p.o. Ifforts to reproduce results have produced variations 


as high ae 4 divisions of the MPO, with an average error of 1 to 1 1/2 





divisions. This represents an error, or null voltage, of 4 to 6 mv in matching 
voltages. 
A change in the instrumentation was effected just prior to the final 


test flight (See Flight Record - Appendix D). The gyro container was 





secured rigidly to the aircraft structure during flight in order to eliminate 
recording vibrations setty by the gyro container mount and the Lord shock 
mounte. This modification resulted in removing an undesirable 10 cps 
transient from the records which occurred whenever the aircraft was respond- 
ing to input pulses. This change made it possible to make the shock mounts 
ineffective during flight, and yet to free them when desired in order to 


protect the gyros from landing ehock. 
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ADJUSTABLE MOUNTING 
FOR GYRO CONTAINER 
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APPENDIX ¢ 
ANALYSIS OF A DIVERGENT TRANSIENT RESPONSE 
References (2) and (3) point out that the approximation of the Fourier 
Transform is applicable only to stable systems. This can be easily seen fron 
a consideration of the transient response of a system described by a 


characteristic equation having a positive reoot. For such a response, the sum 


- jndrw 
r q(t. 2 
nem Vee... — > 


é 


will not converge. This condition is found in the lateral motion of many 
airplanes, where the positive root is associated with the motion called 
epiral divergence, 

Since the B-26J used in the present investigation did not exhibit this 
type of response, a simple dynamic system illuatrating divergent motion is 
analyzed © wus appendix, using the finite pulee method. 

Consider the system schematically pictured in Fig. C-1. The mase, m, 
in this figure is aseumed to be supported in the gravity field by a rigid, 
masslees,rod. The Laplace Traneform of the characteristic equation for the 
emall angles of motion coneidered in Fig. C-1, (Ref. 8), is 

oats 6 2.6 ime) © = 0 (C~1) 
Assume that a forcing function is applied consisting of a unit triangle 


of force, UT(t - At) which is described by 


F= 0 0 >" 
1 < Z 
Ar 
Fu 2- 8 t Ar £ ¢ £& 2d 
Ar 
FF =O t > eat 


Ar = 0.1 second 











Oe wi oe 
- a 


> -@* Ses 
[mL +o a - meloek (Lt)[ ( 


ie _, 7 





(LT)(UT(t - Ar)} = 










Solving for @ (s), 
_ .-O78,2 
© (s) = 5 > ¢ 
(Ar 9 )(ml )(s t as y = 5) 






In this investigation of the B-25J lateral motion, the response 


was measured was the angular rate in roll or yaw. For this reason 





time rate of change of the angular position ® will be considered 


in this example. The rate of change of this output, » is given by 


: _ ate 2 
@(s) = hae 
Ar ol s(s tea- 





Factoring the characteristic equation, 


2 
ae ee ee 
ap-\ 2 J 


Let 
°° ez 
ace > 
‘ | 4m q 
yy) 7 
ae ons 
= 4m 
Note that a is diways positive and that a and b are positive, real 
numbers. 


By appropriate substitution of a and b, the equation of motion, eq. (C-4) 


can be written, 


-A 2 
O (ya: UE (0-5) 


Ov ml s (s - ad(s + d) 
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The performance function (PF) » @ } is defined ae 


Lt) @ 
(Try 67" i rit 


Then, 


= 


; 1 
FFiy, 0)" nl Go- alle od) ae 


A value of 0.1 was assigned to a, and a walue of 2.0 to bd. The value 
of 0.1 for a ia such that the amplitude of the divergence will double in 6.93 
seconds. Some airplanes have this mich divergence at high angles of attack. 
The value of 2.0 for »b is somewhat small compared with that for the highly 
damped motion of most airplanes. The major consideration in the choice of 
these roote was to choose values that would reduce the oalculations as such 
aes practicable, and yet that would retain some semblance to those encountered 
in airplane lateral performance functions. 

These choices were consistent with the choice of mand / as 1 and 2 
reepectively. The units are those of any consistent set of physical unites. 

Subatituting these numbers and letting 

s = jo 


the performance function becomes 
; i 
Fry @) 7" 2 Tw 0.1)(jo + 2) (c-7) 


This performance function is plotted in Fig. C-5, 

Returning to eq. (C-5), by the inverse Laplace transformation the time 
response was obtained and plotted in Fig. C-2. 

In Fig. C-3, the portion of the curve of Fig. C-2 beyond 5 seconds was 
plotted on semi-log paper, thus enabling determination of the divergent root. 
The straight line was extended to zero time and the intercept found. Having 


thus found an expression for the divergent component of the response, it was 








subtracted from the total response leaving the remainder shown in Fig. C-4. 
The integral of thirx remainder is convergent and its approximate Fourier 


Transform is 


(FT) @ (y(t) = (Fryor(s}  z @(t_)e nam 


on) Mes lee... 
Since the system is linear, the performance function determined from 


the response curve can be expressed as the sum of components, i.e., 


(LT)(q (t)] * t) (FT)(q_ (t)] 
yo) 7 * one 
(ryt (RT)Ca (td) on) | 
_ ara allied 


The last two terms on the right side of eq. (C-8) were determined by 
approximating the Fourier Transform by triangles. The first term can be 


determined since the analytical expression for a, (t) is known, and 


divergent 


since 


(LT)Eq, (t)) 7 (Ft)Cq, (t)] 


Thie analytical expreseion for the @ivergent component of the response as 


determined from the time response is 


Gra (tay] = ew 


a — a 
[F, @ J(dtv) (yy ~ 9.1) (Fr)LUT(t) Je O™ 


From this, the amplitude ratio and phase angle can be calculated for any 
given ® , 

The desired performance ‘function can be obtained by performing the 
addition of complex quantities indicated in eq. (C-8). 

Figure C-5 shows the comparison of the original performance function 


with the performance function obtained by analysis of the time response. The 
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FIG. C-1. SCHEMATIC DIAGRAM OF SIMPLE DYNAMICAL SYSTEM 
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APPENDIX D 
DATA REDUCTION AND EVALUATION 

Te method used to obtain a desired performance function, knowing the 
time response of a system to a known input pulse, is outlined in Chapter l. 
In that chapter the method of approximating the mathematical expressions 
graphically (within the desired degree of accuracy) is covered. Reference 3 
shows the method applied to a critically damped second order systen, and out 
lines a method of attack to be used when more lightly damped oscillations 
are present in the time response. The process as above outlined has been 
applied to the oscillograph test records obtained in this investigation. 


Te oscillograph records provided the following information: 


Input Pulses Output Response 
Aileron deflection, a (+ defl) Angular vel. in yaw and roll 
Aijleron deflection, Pia (— defl) Angular vel. in yaw and roll 
Rudder deflection, 6. (+ defl) Angular vel. in yaw and roll 
Rudder deflection, oF (- defi) Angular vel. in yaw and roll 


Sign conventions used in designating pulses correspond to NACA usage, i a 
(+) Aileron defl. - right aileron depressed 
(+) Rudder defl. - rudder trailirg edge moved to the left. 

- Reduction of the above listed recorded data would yield the following 


aircraft performance functions: 


1) (YF (2) (PF 
CO) CPA) (8 ow) Jeay(a., Wy) 


(3) (PF (4) (PF 
CPP Cay (a, Wy) Vaya, W,) 





The time available for reduction of data was not sufficient to allow 
all these performance functione to be determined. Therefore, since this 
project is primarily concerned with the application of the pulse technique 
as a means Of determining the performance function of any dynamical systen, 
and since knowledge of the actual dynamics of the test camber is of 
secondary importance, the records were analyzed to achieve the following 
objectives: 

(1) To determine the degree of accuracy or uncertainty in repeating 
results. 

(2) To compare the results of determining the performance function by 
slightly variant graphical means. 

(3) fo compare the performance function as obtained by the pulse 
technique with the calculated performance function as determined from the 
aerodynamic derivatives obtained from wind tunnel teste. 

(4) To reveal any inherent uncertainties or limitations in using the 
pulse technique to determine performance functions. 

To obtain the data necessary to investigate the above objectives, the 


following analyses were actually made: 











TABLE 1 
Direction of 
Run No. (PFI r } Method Pulse Input 
4665 (PF]+, 40] . Graphical (+) defl. 
a X 
4661 (PF Ir, Ww) Graphical (~) defl. 
a 
4662 (PP]rg w_ Graphical (-) defl. 
a . 
4665 (Pre v,] Graphical (+) defl. 
& 
4662 (FF] "i Graphical (+) defl. 
. a 2 
4665 (FF), Ww.) Numerical (+) defl. 
a 2 





In determining the performance function graphically, three different 
methods were utilized to determine the amplitude ratio and phase angle at 
each selected forcing frequency. The first method utilized the general 
method as outlined in Ref. 3. that is, the graphical summation of input and 


output vectors as indicated in the following expression: 


- jnw Ar 
(FF) = d 
~Jnw,dt, (D-1) 


2 Van) 't,e (FT}(UT(t)], 


Hn ww 152.0484 050 ee 


Figure D-l shows a sample plot of this method for w, = 1.25 rad/sec using 


f 
run number 4665, for determining the performance function CPF] ny ce Ww) In 
this plot, the effect of the lightly damped sinusoid on the phase cain and 
amplitude ratio is apparent. Hach cusp represents one half cycle of 
oscilletion. In order for the amplitude to be attenuated to 0.05 of the 
original amplitude, the plot must be continued until a total of approximstely 
6.5 cusps have been plotted. The majority of the oscillograph records are 
not long enough to allow more than three or four cusps to be plotted, sae is it 
practicable to obtain transient responses of this length with a usable 
uncertainty level. Thus the AR or PA, or both AR and PA are subject to 
error due to deleting the vector summation representing several additional 
half cycles of the oscillatory mode. Furthermore, the desired matching of 
the curve required that the response be divided into 0.1 sec increments, 
Pei r ing that as high as 84 vectors be plotted (8.4 sec) at each forcing 
frequency —- a tedious job resulting in a rather high degree of uncertainty. 

In order to simplify the graphical method by requiring fewer vectors 
to be plotted to determine the AR and PA at each frequency, and to obviate 
the error incurred by omitting two or three half cycles of the oscillatory 


mode, a second graphical plotting system was used. In this system vectors 


were plotted in identically the same manner as previously discussed until 
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the time responee was reduced to a remainder which was solely (for practical 
purposes) a pure damped sinusoid. At this point a single vector was added 
to the response vector summation to represent the contribution of the long 
"*tail’’ of the sinusoid. 

Figure D-2 is a sample plot of this method at the same forcing frequency 
and run (No. 4665) as used in Fig. D-1. The accuracy of this method depends 
on the accuracy with which the characteristics of the damped sinusoid may 
be determined. The characteristics were determined by the method demonstrated 
in Fig. 8 of Hef. 3. Knowing that the oscillatory motions in both yawing and 
rolling angular velocity have the same frequency and damping ratio 
(see eqs. (D-10) and (D-11) several records of both motions were used to 
measure the period and peak amplitude ratios, and average values for the 
parameters, DR and Wa’ were thus obtained. The logarithmic decrement curves 
shown in Fig. 7 of Ref. 3 were used in determining DR. The records agreed 
quite closely in the period of the oscillation, but determination of DR 
required the average of several records, giving the following characteristics: 

DR = 0,15 
W = 1.45 rads/sec 

In the case of the rolling angular velocity response, it was found 
that the response was nearly a pure sinusoid after the third peak overshoot. 
(See Fig. D-5). As a consequence, the general performance function equation 


was rewritten, designating the time at the second peak overshoot as oe and 


-jnw_ Ar 
the ordinate at . ar 1° ( (x Vout) Tn)? ¢ 0 crmoe(0)),1 1.2 


ao it ayn ane ; Jnw 


((z Xan) (Tq)® ony crANLONe I 


1,< ene 


~ jnw At ) 
( (z q, (7 )e a & O° praycor(s)],) a a 
~Jjnw At, 1 
(x (7) ) (FrrjCur(t) 
ee yh “7 4 = 1,2,... (oe 





(zero time lag, Tt)? ie assumed in the preceding equation) 
The second term in the numerator may now be replaced by a single vector 


representing the contribution of the damped sinusoid from a to infinity, 








- jnw Ar 
yielding, Facute * ° Crritur(t)), 
| ea ge a Ee... 8 
LPP] (gS ley ~Jnw At, : 
x Van) (Te (FT)(UT(t)], 
n= Me ig tene ie 
~Spw Ar 
({FT]lq autre ee . ie (D-3) 
-Jjnw AT 
2 Vann 6 e ; (Frj(ur(t)), 23 
n = ie eVMnece 





where Or and 4r, are the time intervals used in the output and input 


4 
respectively. 

Any time later than the time at which the response becomes a pure sinusoid 
may be selected as the point at which the remaining oscillatory motion may be 
represented by the Fourier Transform of a damped sine or cosine wave. However, 
it i@ advantageous to use a *’peak point’” or a point where the sinusoid 
crosses the axis in order that the phase angle in the sinusoid expression may 
be eliminated. In these two cases the term is reduced to the following Fourier 


Transform expressions: 


t at peak) 
(t,, p 











me ~Jpw Ar 
CFT] (osc) = L¥TJq., 6 cos wt | ° 
= = ~ Jpw Ar 
1 1 é 
q , (DR) 2 * By’ tan (FR) - tan 2(DR) (FR) e 
“|W 
° [1 - (FR)°}° + [e(pR) (¥R)]° (D4) 
(t at axis) 
. ~(DR) (wt ee, 
[FT )a0 (onc) = [FT Ja, a gin wit e 
to © ~Jpw Ar, 
: , 5 6 (D—5) 
n 
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The time, Ms should be chosen to fall on one of the ordinates determined by 
the choice of Ot to avoid disturbing the area approximation originally 
established by the isosceles triangles. If a peak maximum or minimun is 
selected as ty the area approximation by triangles is interrupted even 
through SS may lie on one of the ordinates of the fitted triangle. In this 
case, a right triangle term mest be added to the Fourier Transform of the 
output if the resulting output vector is to be exactly correct. However, 
the initial amplitude of the cosine wave starting at Me was small enough 

in the runs analyzed in this project to render the right triangle term 
ineffectual in changing the AR or PA of the performance function. The 
method is demonstrated in Fig. D=4. 

Figures D-1, D-<2 and D-3 compare the two graphical methods used in 
determining the performance function, LPPJCay(s., ¥,)” Figure D-l shows 
the AR and PA determined by plotting the maximum number of vectora 
available from the oscillograph record. Figure D-2 shows the same AR 
and PA at the same forcing frequency as determined by representing the 
output divided in two parts in the method just described. Figure D-3 
shows a comparison of the AR and PA resulte determined by each method. 

It was found that the maximum error in the first method, Fig. D-l, appeared 
at the resonant frequency. 

The third method used for determination of the performance function 
from pulse response data is the method described in Ref. 3, Section IV. 

The oscillatory mode was determined ag previously described, and 

eubtracted from the output response from t = 0 to ©. As shown in 

Fig. 9 of Ref. 3, the remainder of the output response is no longer 

zero at t = ae and a term must be added to account for the right triangular 


area at the time origin, which is not included in the area approximated by 


=-9R=— 





the isosceles triangles approximating the remainder of the response. The 


resulting expression is that giwen by eq. 40 of Ref. 3. 


_ FN, (ope) 1) * Ag (rem) (Tq) MPEIC CURE) CH) 


~Jnw Ar, 
; : bo, .Mta)n?® (FT}{ (UT) (t)), 


eer May. dea! 


- jnw OT 
y(r )Je fi" (FT]((UT)(+)}, 


- jnw At 
r (a, (7 )Je ons [FT)((UT)(t)), 


tiq 


o(rem =: 1.2.08 


N= 1,2, 00. eee 

Figures 3-1 and 32 show the amplitude ratio and phase angle 
respectively of the performance function for aileron deflection input and 
roll response output. The phase angles shown have been corrected for gyro 
phase lag as obtained from the Fig. D-13. Undamped natural frequencies and 
damping ratio for both gyros are found in Appendix B. In figures 21 and 
3-2 three experimental performance functions are shown which are the result 
of analyzing graphically three different oscillograph records (rune 4665, 
4661 and 4662). The uncertainty in reproducing results may be observed 
from these plote. The curve shown represents the average value of the 
three experimental points at each frequency. The greatest deviation from 
the average occurs at We = 1.25 where the maximum wariation in AR is 5.9 
percent. The average variation is 1.6 percent. As noted in the sample 
oscillograph records, Figs. D-5 and D-6, no two input pulses were identical. 
Consequently, each hep oni is different which necessitated independent 
analysis of each record. The reliability of the pulse input technique in 
reproducing results appears to be excellent, judging from the results of 
this limited investigation. 

From the experience gained in use of the analyzing systems previously 
outlined, the following advantages and disadvantages were noted in applying 
these systems to time records having a long ’’tail’’ due to a lightly damped 


sinusoid. 
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The only remaining curve shown in Figs. $-1 and 22 which has not been 
discussed is the calculated performance function using wind tunnel stability 
derivatives for the B-25. Due to differences between the model used in the 
wind tunnel tests and the actual airplane, these derivatives have been 
modified to more closely represent the true derivatives of the aircraft. 
The derivatives that were not available from wind tunnel teats were estimated 
on the basis of previous engineering experience with comparable aircraft 
configurations. For use as a means of comparison in this project, the 
derivatives were corrected to include the values of moment of inertia about 
the X and Z axes as found for the B-25 used in the testing configuration 
(See Appendix A). 

The derivatives as used to determine a calculated performance function 


are as follows: 











ac 
L tle Cau = -2.71 1/sec 
P o(Be) Tix 
20 
oc 
L = n= 2 = +,673 1/sec 
TH 
0G Py 
L i "oh = - 00964 1/ft-sec 
ac 2 2 
a Sw eee. 2 
lL, * 56 2, 3.175 1/sec 
0c 
2 2 LOO oe a3 l/pae’ 
Sa 3 be 21, 
dc 
f° SoU 
Ns a = -.10 1/sec 
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0(5e) ZZ 
Qo. ss 
N ee = +,00648 1/ft-sec 
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do 


/?USd 
. 2a Z 
KF = 7iL> * ~ 48 1/sec 
o (SR) ZZ 
oc 
/? 8U 
as = a) 2 [= = -.13 1/sec 


7 oo ? spu’ 


2 
-1.66 1/sec 
é6r 3 * 2100 


Alt. = 10,000 ft 


Air epd = 175 mph 


2 
Iyo = -1930 slug ft 
2 
“yy = 63,000 slug ft 
I = 120,000 elug et” 
ZZ 
Weight = 26,000 lb 


The sign convention for the derivatives and equations of motion and 
the symbols representing displacements, moments, velocities, etc., are in 
accordance with NACA designation system. Positive directions along the 
axes, viewed from the origin, are forward, right wing, and down for the 
X, Y, and Z axes respectively. 

The following assumptions were made in writing the equation of motion: 

1. Damping momenta, dire to control surface motion are negligible. 


2. The wind axee and the aircraft reference axes are considered to 
be superimposed. 


3. No gyroscopic effects were considered. 
The equations of motion are; 


a. Summation of forces along Y axis 


om 


é 
m + mr = rte ase” (D~7) 
bd. Summation of rolling moments 
- OL, , bi, ob 4, Oh 
Tre PiT Oy * Op?” Or * Mr "Gee an (D-8) 
c. Summation of yawing moments 
- OL, Ly ok a: 
Ms? < 5e* op PO toor® tae? * eS oe 
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Dividing each equation by the appropriate moment of inertia or the 
mass, taking the Laplace transform and solving for p and r, the following 


performance functions were obdtained: 


-3.175 8 (s° + ,0868 + 1.6) 


PY ‘it 
PING My)” (52 4 oar7 s+ 1.78)(8 + 2.9)(8— .0013) (D210) 


where p = Wy. 


1184 (§ + 3.51)(8 + 1.097)(S - .837) __ 
LPFJ(g w,)* 


(3° + .377 S + 1.78)(S + 2.9)(8 - .0013) (B11) 
where r 2 Woe 
The calculated performance function LPF (g ¥, ) shown in Figs. 3X1 
and 3-2 have the same pared shape as the shane performance 
function. The misalignment of the resonant frequencies of the two plots 
ie indicated in the difference in characteristics of the oscillatory modes 
of the calculated and experimental performance functions. The characteristics 


for the calculated and experimental oscillatory modes are: 


DE W 
n 
Experimental  .15 1.48 
Calculated 0141 1.36 


Since the oscillatory mode is the predominant mode in establishing 
the peak values of the performance function, the difference in the 
magnitude of the two peake should represent approximately the sane 
percentage difference as indicated by the difference in their damping 
ratios. From Fig. &1 the peak of the experimental performance function is 
93.8 percent of the peak value of the calculated performance function, while 
the corresponding ratios of the two DR is 93 percent. 

The calculated performance function shows a very small divergent root. 
However, it is difficult to state whether the B-25J actually is spirally 


divergent or not. The oscillograph time records do not indicate a divergence. 
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However, in order to determine whether such divergence exists, it would be 
necessary to take a record of at least 15 seconde to determine the presence 
of this unstable mode. During this period small perturbations in roll or yaw 
due to disturbances of the air mase would probably nullify the accuracy of 
any such record. In any case, for practical considerations, a root as small 
as the calculated performance function indicates, whether positive or 
negative, would have negligible effect on the performance function. 

The sensitivities were compared using the following assumptions: 

1. The unstable root encountered in the characteristic equation of 

the calculated performance function was considered negligible 


(0.0013 = 0), 


2. The experimental performance function contains no divergent roots. 
(An assumption implicit in positioning the base line.) 


Using the first assumption, the sensitivity for the calculated response was 
determined by rewriting the performance function as a sensitivity mitiplied 
by a frequency function. The same assumption applied to the experimental 
curve gives the ratio of the area under the responee curve in the time 
domain to the area under the input curve in the time domain as the sensitivity 
of the experimental performance function. 

To determine the area contributed by the oscillatory mode of the 
response extending beyond the limits of the oscillograph record, the area 
was measured out to t = t (t,, again representing the point at which the 
response becomes a pure damped sinusoid) and the curve beyond time ms was 
integrated analytically and added to the area from t = 0 to : The 


expression thus becomes t = oo 


~DEW t 
[ZL q (+_) dr_] + fq (r) e coswt 
60 jn ae : Pe 6 
2 
q, si At 


S(ay(6 Ww 
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E’ oxperidsiital i 


n = Ec ears 
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The sensitivities thus determined are: 


S -.992 1 
(A) (8, "x (cate) bees 


S -.945 l/sec (Run Wo. 4666 Input = + 6 ) 


-1.080 1/sec (Run No. 4665 Input = + 8) 


-.985 1/sec (Bun No. 4661 Input = - 8) 
==990 1l/sec (Run No. 4€62 Input = - 6) 


S =-1.0 1/sec 
(A) (6. Wx) (exper) (av) 


An additional sensitivity for comparison with the above may be obtained from 
flight teet data for the B-25, utilizing curves of 32 vB . (Ref. 4). From 


these curwes, the sensitivity is, 


g =-90 1 
(A) (8 Wy) (F) (test data) ies 


The comparison of experimental and calculated sensitivities are in 
close agreement, as is the check with flight test roll rate data (Ref. 4). 

The performance function for yawing velocity response to aileron input 
is shown in Figs. 23 and 34. The vawing velocity response as shown in 
Fig. Ds appears to be an excellent example to illustrate the limitations 
of the pulse input technique. As noted in Fig. De-5, the response is 
essentially oscillatory, the first lobe (positive) indicating the adverse 
yaw effect in which side slip has taken place to give the high righting 
moment responsible for the second lobe (negative) being of greater amplitude 
‘than the first lobe. The pattern of this response is noted in detail since 
its effect in the graphical analysis is important. 

The type of output response pattern found in yawing velocity for aileron 
input shows the same summation pattern at each of the low frequencies. Com 


paring Fig. De? with Fig. De2, it is seen that the end point of each 


els 












ummation is well away from the origin, and that small errors made in 
plotting will offer very little error to either the AR or the PA of the 
total output vector. Also, the comparison at these low frequencies with 

the AR and PA of the calculated performance function is satisfactory both 

in yaw and in roll. (See Fig. 3-1 and 32 for roll, and Fig. 3-3 and 3-4 
for yaw.) However, as the frequency increases, and WAT increases corres- 
pondingly, the vector summation of the yaw vectors begin to ‘'’wind up,’’ 

the end points never getting very far from the origin. Note the difference 
in the summation as shown for W, = 1.5 rad/sec in Fig. D-? and the summation 


as shown for W, = 5.0 rad/sec in Fig. D-8. This clearly indicates that any 


fT 
uncertainty entering the summation, regardless of its nature, will impart a 
mich higher uncertainty effect to the PA of the output vector (end-point to 
origin) than in the case where the end-point is well removed from the origin. 
The uncertainty of the results thus are increased by two ere: as the 
frequency increases, 

1. The large relative size of the last vector in the summation in 
comparison with the resultant output vector, and 

2. The accumulated errors from measuring ordinates from the record. 

Consequently, each system of analysis as covered in Table 2 meat be 
analyzed and tested to determine the relative uncertainty in approximating 
the Fourier Transform under conditions where the response is composed in the 
main of an oscillatory motion. The general effort is to establish the system 
of analysis that will produce the fewest inaccuracies in approximation and the 
fewest uncertainties in plotting. 

Inveatigation of the method listed in Table 2 as (2), the one used 
most extensively in determining the response vectors for the performance 
function in roll for an aileron input, certain inherent uncertainties are 


noted. Again, it is emphasized that these uncertainties become increasingly 
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important as the response becomes more predominately oscillatory. These 
uncertainties affecting the magnitude and phase angle of the output vector 
are listed in order of the importance of the error they present. 

1. Uncertainty in determining the time at which a mode has reached its 
peak value on the time response, or the time at which the output time response 
crosses the base line. That is, correct determination of the location ef oe 
Figures D-S and D-6, actual oscillograph records, show the cause for this 
uncertainty. The records are of too short duration to show more than a maxi- 
mum of one full cycle of the pure damped sine wave. Therefore, although the 
DR and a of the oscillating mode may be well established by comparison of 
many records, the low flat=-topped lobes and the low slope of the curve, plus 
the fact that the noise to signal level is becoming quite high at low 
amplitudes, all tend to make the selection of the exact point at which the 
curve peaks or crosses the base line somewhat uncertain within about +0.1 
second. The effect of such an error can be observed best at high frequency. 
If We = 6 rads/sec anddt = 0.1, and the final output vector 1s approximately 
the length of the vector representing the oscillatory mode, an error as high 
as 30 degrees in PA may exist. As the end point moves further from the 
origin, this possible maximum error reduces accordingly. 

2. Uncertainty in determining DR and Wo of the oscillatory mode. 

3. Uncertainty in determining the amplitude of the oscillatory mode. 

On the basie of the above uncertainty at high frequency in adding a 
vector representing the Fourier Transform of the oscillatory mode, it may 
be concluded that the greater the time, oa the less the uncertainty. This 
may be observed in eq. (D-3). To simplify that expression consider 


Ar = At., then 
0 i 


= Ope 
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E = lic, 
LPP] (q,a.) ~ ~jnw At, 
od (T Je Wetec po is oa 
~JPAT We 
| [orn] 
~Jaow Ar, 
C(rrj(ur(t)), z a, (7 Je mn alee ue 


The larger a becomes, the greater the percentage of the graphical analysis 
that is accomplished with the same degree of approximation in both numerator 
and denominator. If the oscillograph record were long enough, the best fesults 
at high frequency would be obtained by continuing summation of the output 
vectors until such time as the end point is clearly indicated. This allows 
the same degree of approximation in obtaining both the output vector and 
the invut vector throughout the entire summation. However, in addition to 
the tedium of plotting a large number of vectors, uncertainties are intro- 
duced by the oscillograph record in picking off ordinates well out in the 
time respcnee. The noise to signal ratio becomes high, and drift of the 
baseline appears to be one of the evils of the long record. 

The most uncertain of the graphical methods for use with the highly 
oscillating type response is the method which subtracts the pure 
oscillatory mode from the entire response. In this case, the subtracted 
oscillatory mode represents a high percentage of the total response. This 
large portion of the response is represented by a single vector which is an 
exact Fourier Transform. Consequently, the output, in the extreme case where 
the response is almost entirely oscillstory, is represented exactly, and the 
input is approximated. The result 18 that at high frequency, the output 


reaches a certain fixed phase angle, due to the predominance of the oscillatory 
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mode, and the input continues to shift phase angle rapidly. Therefore, 
it appears that at very high frequency the oscillatory part of the output 
will eventually lead the input. The lagging phase of the remainder should 
theoretically compeneate in such a manner that the vector representing the 
total performance function will have the proper phase angle. But, the 
vector representing the remainder is very eamall, and therefore, although 
its phase lag may be increasing in a like manner to the denominator (input), 
its contribution to the phase angle of the vector summation is very small. 
Furthermore, the additional manipulation of the time response data in 
obtaining a breakdown into components offers added uncertainty in the vector 
summation. 

Sample records are provided to show the vector summation using method 
(3) in which the oscillatory mode is removed from t = 0 to t = °© and 
represented by an exact Fourier Transform, and method (2) in which the 
Fourier Transform of a sine wave starting at ‘ is added to the output 
summation for t = 0 to t = “i Figures D-9 and D-10 respectively show these 
two cases. Note that in Fig. DB the vector representing the oscillatory 
mode is the largest of the three vectors comprising the total response 
vector. Furthermore, it has attained its maximum lagging phase angle to 
within about 5 degrees. The phase angle from the performance function at 
We = 6 by this method is -49 degrees (corrected for gyro phase angle 
from Fig. D-14). Figure D-10 shows the resulte of using the same approxi- 
mation in both input and output over the greater part of the response, method 
(2). The addition of the vector representing the sine wave summation is 
nearly negligible in its effect on the total response vector. The phase 
angle at this same frequency by this method, (2), is -87.5 degrees. 


As previously discussed, adding the eine term involves uncertainty in 


properly locating ve This error may be reduced to a negligible size or 
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eliminated by the following means. Select es where the output crosses the 
base line and observe the magnitude of the vector resulting in comparison 
with the fan! tithe of the last vector in the summation from t = 0 to t = Ege 
If the vectors are approximately the same size, no appreciable error will 
result from an error in locating eas If the vector representing the added 
eine term is large compared to the last vector of the summation, - should 
be moved to the next null point of the time response, which will decrease 
-W (DR) (T/2) 
the magnitude of the sine wave amplitude by the factor e ; é 

It is the opinion of the authors that the added sine wave method of 
component summation is subject to the least inaccuracies, of the methods 
investigated, in determining the performance function of a response which is 
essentially oscillatory. The method affords the same degree of approximation 
in both the input and response, until such time as the response is subjected 
to inaccuracies of noise and instrumentation. At this point (t the small 
remainder is represented by an exact transforn. 

Referring again to Figs. 3-3 and 3-4, the effects of various analysis 
methods in obtaining the performance function are shown. The method numbers 
are those listed in Table 2. As previously noted, the variation in PA from 
the calculated performance function at high frequencies is very high, and it 
is difficult to explain the variation in terms of the dynamics of the air 
plane. The assumptions listed in setting up the theoretical equations of 
motion for the test vehicle, limit the phase shift to 90 degrees of lag at 
high frequencies. However, refinements which were excluded by these 
assumptions may add about 10 to 15 percent additional phase lag at high 
frequency, as shown in the theoretical performance function found in Ref. 7. 
The performance function applies to the B-26, an aircraft similar in dimensions 
and mass configuration to the B-25J. The experimental PA, at high frequencies 
are in disagreement with the calculated performance functions determined both 


for the B-25 and the B-26. However, the very low AR at the frequencies of 
Sie — 





uncertain phase angles makes PA uncertainty relatively unimportant in 
determining the performance functions. 

The scattering of points in the doubtful phase angle range shows the 
effect of varioue methods of determining the transform of the response. It 
was found by halving the Ors the phase angle was negligibly affected. 
Compare Fig. D-8 with Fig. D-ll. This indicates that W Ar has not reached 
a critical value. (Ref.3). 

The sensitivities for this performance function were determined in the 
same manner as used in the performance function for rolling velocity response 
to input aileron reeponses. 

Ae previously noted, the calculated performance function, (eq. D-1l), 
has a very small root in the denominator. Such a small root either positive 
of negative would not be apparent in the time response and would be implicitly 
considered to be zero by the action of drawing in the base line. To be 
consistent for comparative purposes, the same assumption mat be made with 
the calculated performance function. 

Thue, to get a valid comparison let jw - 0.0013 = jw in the calculated 
performance function. Then by writing the performance function as a 
sensitivity multiplied by a frequency function, a sensitivity that can 
ve compared with the ratio of the areas under the experimental output and 
input curves is obtained. The comparison is: 


alls, %) = 0.209 1/sec 


cale 
a ‘ + 6« f t 
SAEs 4) 167 1/eec (from run no. 4665 , inpu ) 
a Z-exper 
= .17 ir . 4662 —6 { t 
Sats wv 1 172 1/sec (from run no , inpu ) 
a Z-exper 


In general, the results obtained bring out several points concerning 


the application of the pulse technique as a means of determining the 
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performance functions of dynamical systems. Results obtained in rolling 
velocity response for aileron inputs showed very close agreement with the 
theoretical performance function throughout the entire range of frequencies 
investigated. The comparison of time response patterns or conformations 
indicate the resvonse pattern involving the least inaccuracies in analysis 

is one which responds quickly, reaching its maximum amplitude in the early 
part of the time response. Under these conditions, the end point will usually 
ve well removed from the origin, and uncertainties in added vectors represent~ 
ing the Fourier Transforms of oscillatory modes, right triangles, etc. will 
have a minimum effect on the AR and PA of the performance mnct low until 

very high frequencies and small AR prevail. 

A potential source of error lies in establishing the true baseline for 
each response and input time record. An error in the location of this line © 
produces a double error in the results in that as much false area is added 
to one side of the line as is subtracted from the other side of the line. In 
the case of records having the long tail characteristic of a lightly damped 
oscillatory mode, the envelope should be drawn in determining the base line, 
particularly if there is any drift in the center of the response record. This 
is demonstrated in Fig. 8, Ref. 3. In any case, measurements from the refersnce 
trace to the center-line should be made at short intervals and connected, 
rather than passing a line between two measured points at the extremities of 
the record. This corrects for warping of the record incurred in the drying 
process. If nominal care is exercised in laying out the baseline, and in 
picking off ordinates from the record, the uncertainties in reproducing 


results will be satisfactorily low. 
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FLIGHT TEST RECORDS 
FLIGHT NO, 1 
Date: 28 July 1949 Duration: 2.5 hres 
Airplane: B-258J No. 44-30328 
Crew: 
Pilot: C.0O. Bostrom 
Co~pilot: F.H. Michaelis 


Observers: J.C. Wootton 
A. A. Hollander 


Purpose of Flight: 
Instrument shakedown 
Airplane configuration: 
Gross Weight - 26,100 
C. @, Percent MAC - 27.4 
Instrumentation installed for recording elevator position, rudder position, 
rate of roll and rate of yaw. 
Flight Procedure: 
(1) Trimmed aircraft at 10,000 ft PA and 175 mph IAS. 
(2) Observed action of oscillograph of input and response traces 
using test pulses of positive and Pot tic eth ileren and rudder displacemente. 
(3) Obtained records of positive and negative ailerons and 
rudder pulses. 
Results and Comments: 
Frequency regulation of the inverter was very poor. Frequency varied 
from 380 cps to 405 cps. Noise level on the oscillograph records was s80 


high that the records were not satisfactory for reduction. 
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FLIGHT NO, 2 
Date: c August 1949 Duration: 1.5 hrs 
Airplane: B-25J No. 44-30328 
Crew: 
Pilot: C. Collins 
Compilot: J.B. Bain 


Observers: F.H. Michaelis 
J.C. Wootton 


Purpose of Flight: 

To obtain oscillograph records of airplane response to aileron and 
rudder inputa. 
Airplane configuration: 

Gross welght.- 26,100 

C. G@. Percent MAC — 27.4 
Yiight Procedure: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Ovserved action of oscillograph traces for teat pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 


No records obtained due to malfunctioning of recording equipment. 
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FLIGHT NO, 3 
Date: c¢ August 1949 Duration: 1.5 hrs 
Airplane: B-2&5J3 No. 44-30328 
Crew: 
Pilot: C. Collins 
Compilot: J.B. Bain 


Observers: F.H. Michaelis 
J.C. Wootton 


Purpose of Flight: 

To obtain oscillograph records of aircraft response to aileron and 
rudder inputs. 
Airplane configuration: 

Gross weight - 26,100 

C.G. Percent MAC - 27.4 
Procedure this Flight: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Observed action of oscillograph traces for test pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

The records obtained had an unacceptable noise level, but were somewhat 
improved over preceding flights. Input pulses, while acceptable, were too 


long for ease of data reduction. 
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FLIGHT WO, 4 
Date: 5 August 1949 Duration: 1.50 hr 
Airplane: B-25J No. 44-30328 
Crew: 
Pilot: Capt. Warwick 
Co-pilot: J.B. Bain 


Observers: J.C. Wootton 
F.H. Michaelis 


Purpose of Flight: 

To obtain oscillograph records of airplane response to aileron and 
rudder displacements. 
Airplane configuration: 

Gross weight - 26,100 

C.G. Percent MAC - 27.4 
Replaced faulty inverter, and installed new Lord shock mounts on gyro can. 
Flight Procedure: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Ovserved action of oscillograph traces for test pulses of 
positive and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulees. 
Results and Commente: 


Flight discontinued due to wiring difficulties. 
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FLIGHT NO, 5 


Date: 5 August 1949 Duration: 


Airplane: B-25J No. 44-30328 
Crew; 
Pilot: Capt. Warwick 
Compilot: J.B. Bain 


Obserwers: J.C. Wootton 
¥.H. Michaelis 


Airplane configuration: 

Gross weight - 26,100 

C. G. Percent MAC - 27.4 
Changes since last flight: 

Loose connection in 300 volt system found and repaired. 
Flight Procedure: 


(1) Trimmed at 10,000 ft PA and 175 mph IAS. 


(2) Observed action of oscillograph traces for test pulses of 


positive and negative aileron pulses. 


loo mre 


(3) Obtain records of positive and negative aileron and rudder pulses. 


Results and Comments: 


The records obtained were considered to be useable. The noise was of 


insignificant amplitude, and the inout pulses were of the correct size, 


however, the roll gyro picked up a vibration, apparently from its mount, 


signal of about 10 cps which would make data reduction somewhat difficult. 
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FLIGHT NO. 6 
Date: 9 August 1949 Duration: 1.5 Wes 
Airplane: B-25J No. 44-30328 
Crew: 
Pilot: C. Collins 
Co-pilot: J. 8B. Bain 


Observers: F.H. Michaelis 
Pe Suv 


Airplane Configuration: 

Grose weight - 26,100 

C. G@. Percent MAC - 27.4 
Changes since last flight: 

Intalled stop to remove shock mounts while airborne and wedged gyro 
container rigidly to airplane to eldiminate vibration of the mount from 
being recorded. 

Procedure this flight: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS, 

(2) Observed action of oscillograph traces for test pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

The records obtained on this flight were very good. The attenuator 


settings used were: 


Function Attenuator 
Wo S 
Wy 4 
§ 5 
a 
on g 


These attenuations allowed the traces to remain on the paper and gave large 


enough traces to be readable. Eight records were taken and static 
Sale 





calibrations with the master pickoff were made in flight. 
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APPENDIX £ 


GLOSSARY 


The system of notation adopted for use in this thesis has been 


selected becmee of its peculiar usefulness in representing operating 


system performance. 


advantages: 


1. It ia easily learned. 


In particular, it has the following specific 


2. It is adopted to a wide range of situations. 


3. It is built up almost exclusively of characters found on the 


keyvoard of a standard American typewriter. 


4. Any one of the compound symbols of the system is readily inter 


preted without recourse to an extensive glossary. 


A representative list of the primary symbols and those used as sub- 


scripts is given in the following table. 


Primary Symbols 
A 


a ayn Vout? 


(CT) 


(DE) 
(FR) 


(FT) 


XZ 


Angle 


Input - output amplitude ratio fora 
given operating component = Tout)a 


Tina 
Linear viscous damping coefficient 


Characteristic time 


C 


Dasping ratio SS 
2 V ar n 
W 


Frequency ratio = = 
n 


Fourier transform 
Force of gravity/unit of mass 


Moment of inertia of airplane about 
X axis 


Product of inertia 


=132= 





ZZ 


PD May Sou 


t 


im Mas n Sout 


oak May, qeet! 


4(4n) 


Y(out) 
q(t) 


May, Vout) 


~135- 


Moment of inertia of airplane about 
Z axis 


Rolling moment 

Rolling moment due to rolling velocity 
Rolling moment due to yawing velocity 
Rolling moment due to aileron displacement 
Rolling moment due to side slip velocity 
Direct Laplace transform of q(t) = Q(Ss) 
Mass 

Yawing moment 

Yawing moment due to rolling velocity 
Yawing moment due to yawing velocity 
Yawing moment due to side slip velocity 
Yawing moment due to aileron deflection 
Yawing moment due to rudder deflection 


Input - output performance operator for 
a given operating component 


Input - output dynamic phase angle for a 
given operating component 


Input - output performance function for a 
given operating component, where 
Wit + (PA 
st f° an) 


t(in)  Tinlee™ 
Jtw,t i (PA) (out) J 


Vout) ~ T(out)a © 
Input quantity 
Output quantity 


A quantity which is a function of the 
real variable, time 


Right triangle 


Complex variable used in Laplace transform 
theory 


Static sensitivity for a given operating 
component 





At 


(UT) 


X 


Vv 


Time delay 

Finite increment of time 

Airplane forward velocity 

Unit triangle 

Unit right triangle 

Acceleration along Y axis 

Airplane angular velocity about axis ( ) 
Angular forcing frequency 

Natural frequency 

Side force 


Side force due to side slip velocity 


A dot over a variable represente its derivative with respect to time 


(e.g. oO = fo I 


Two dote over a variable represents a second derivative with reepect to 


o°8 2 
time (e.g. ( ) = a () i 
at? 


The following Greek letter symbols are retained because their meanings have 


become thoroughly identified with theese symbols in the aeronautical engineering 


field. 


6 
& 


Total aileron deflection (sum of left and 
right aileron deflection) (positive deflection 
of right aileron is downward) 


Rudder deflection (left rudder deflection is 
positive) 


Angle of bank 
Angle of yaw 


Angular forcing frequency 


-1l|H— 





Modifyi symbols (subscripts 


a Aileron 
(a) Amplitude 
A Airplane 
(div) Divergent 
(in) Input 
(ind sys) Indicating system 
| M Measuring 
0 Out put 
(out ) Out put 
(osc) Oscillating 
p Rate of change of bank angle = Wy 
(rem) Remainder 
r | Rate of change of heading = We 
v Side slip velocity 


=—1o2— 
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